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ABSTRACT 

Future  NASA  plans  to  launch  large  space  structures  solicit  the  need  for  effective  vibration 
control  schemes  which  can  solve  the  uniqtK  problems  associated  with  unwanted  residual 
vibration  in  flexible  spacecraft  In  this  woik,  a  type  of  input  command  shtqting  called 
impulse  shaping  is  examined.  These  shapm  require  only  a  simple  model  of  the  system 
and  are  faMy  insensitive  to  parameter  variations. 

A  theoretical  backg^und  of  inqtulse  shaping  is  presented  along  with  some  insight  into  the 
methods  of  calculating  multiple  mode  sequences.  The  Middedc  Active  Control  Experiment 
(MACE),  an  MIT/NASA  space  shuttle  experiment  is  then  described  as  the  testbed  for  tiie 
experiments.  A  nonlinear  simulation  is  us^  to  verify  the  effectiveness  of  these  shapers  on 
a  model  of  MACE.  Finally,  hardware  experimental  results  are  presented  and  analyz^ 

Some  conclusions  were  derived  from  the  results  of  testing.  First,  it  was  found  that  the 
shaper  was  effective  in  canceling  the  slowest  frequency  of  the  structure,  but  was  less 
successful  in  the  higher  nxxles.  Second,  large  an^e  movements  decreased  tiie  effective 
vibration  cancellation.  The  nonlinearity  involved  with  large  angle  noovements  causes  the 
shaper  to  lose  some  of  its  effectiveness.  Third,  although  methods  exist  which  can  increase 
the  robustness  of  an  impulse  shaper,  ntxie  were  found  to  aid  in  canceling  the  effects  of 
nonlinearities  caused  by  large  an^e  naoves.  Fourth,  the  bandwidth  of  the  controller  was 
found  to  be  linked  to  the  effectiveness  of  the  shaper.  A  higher  barxlwidth  produced  more 
favorable  results. 


Aldiough  many  of  the  inadequacies  of  impulse  sharers  when  applied  to  coix^lex  nonlinear 
systems  are  discussed,  it  remains  a  very  effective  tool  in  ccmtrolling  vibratitm  in  flexible 
space  structures. 


Thesis  Supervisor 


Warren  P.  Seoing 

Professor  of  Mechanical  Engineering 


3 


Aeeession  For 


■TIS  GPA4I 

3^ 

DTIC  TAB 

□ 

Unarmo  ,a:ccd 

□ 

Juatlf IcatioD _ 

By - 

Dl9t,rlbutl''r  / 


Av*:*  1 

I 

Dist 

Special  1 

n 

■gis, 

■  uKoAiiiiL' J 

ACKNOWLEDGMENTS 


I  have  many  things  to  be  thankful  for.  Much  of  what  I  have  received  in  the 
past  two  years  1ms  been  directed  through  the  talents  of  the  people  Fm 
surrounded  by.  This  page  can  only  act  as  a  small  symbol  of  the  gratitude  I 
feel  for  tihe  people  mentioned. 

I  have  to  start  with  my  thesis  advisor,  Warren  Seering.  Thank  you  veiy 
much,  Warren,  for  providing  me  with  Uie  rig^t  balance  of  supervision  and 
autonomy  I  needed  to  finish  my  degree.  You  always  seem  to  know  exactly 
what  to  say  when  I  need  to  hear  it.  Fve  learned  much  firom  you. 

My  education  experience  would  have  been  cheated  had  it  not  been  for  the 
wacl^  group  of  guys  in  my  researdi  group  at  the  AI  lab.  Brian  Avery,  Mike 
Caine,  Andy  Christian,  Jim  Hyde,  Sarath  Krishnaswamy,  Whit  Rappole, 
Lukas  Ruei^er,  Tim  Tuttle,  and  Erik  Vaaler;  thanks.  Thanks  for  putting 
up  with  all  my  stupid  questions  and  thaxiks  for  usually  answering  them. 
You’re  quite  a  bunch. 

I  am  also  in  debt  to  the  MACE  research  group  at  the  Space  Engineering 
Research  Center.  Prof.  Ed  Crawley,  Dr.  Dave  MQler,  Dr.  Mathieu  Mercadal, 
Erik  Saarmaa,  Mark  Campbell,  Dan  Rey,  Carlos  Padilla,  and  Marco 
Quadrelli;  thanks  for  taking  me  in  and  tak^  care  of  me.  You  guys  really 
Imow  what  you’re  doing. 

My  family  has  been  a  constant  source  of  faith  and  support  throughout  my 
years  at  school.  I  wouldn’t  be  near  where  I  am  today  without  them.  Mom, 
Dad,  Bro,  Sis,  Bro-in-Law  and  my  beautiful  baby  niece;  I  hope  you  realize 
how  grateful  I  am  to  all  of  you.  Hey  Eryn,  you’re  in  the  first  chapter! 
Umma,  Appa,  I  love  you  and  this  is  for  you.  I  also  need  to  thank  my  friends 
John,  Jooyee,  and  A^es  for  helping  me  with  the  roui^  transition  to  Boston. 
Thanks  guys,  I  wouldn’t  have  made  it  without  you. 


This  report  describes  research  performed  at  the  Massachusetts  Institute  of 
Technology  Artificial  Intelligence  Laboratoiy  and  at  the  MIT  Space 
Engineering  Research  Center.  Funding  for  this  research  was  provide  in 
part  by  the  National  Aeronautics  and  Space  Administration  under  grant 
#NAGW-1335.  Additional  funding  was  provided  by  the  Office  of  Naval 
Research  under  the  University  Research  Initiative  contract  #N00014-86*K- 
0685. 


4 


TABLE  OF  CONTENTS 

ABSTRACT . 3 

ACKNOWLEDGEMENTS . 4 

CHAPTER  1 

INTRODUCTION . 13 

1.1  Background . 13 

1.2  Input  ^laping . 15 

1.3  Literature  R^ew . 16 

1.4  Overview  of  Thesis . 20 


CHAPTER  2 

THE  TECHNOLOGY:  IMPULSE  SHAPING . 23 

2.1  Introduction . 23 

2.2  Implementing  t^e  Sequences . 24 

2.3  Calculating  iSngle  M^e  Sequmices . 24 

2.4  Calculating  Multiple  Mode  Sequences . .30 

2.4.1  Convolution . 30 

2.4.2  Direct  Solution . 31 

2.4.3  Convolution  vs.  Direct  Solution . 33 

2.5  Finding  the  Most  Efficient  ScKjuence . 37 


CHAPTERS 

THE  TESTBED:  MACE . 41 

3.1  Introduction . 41 

3.2  Project  Motivation  and  Agenda . 41 

3.3  Current  Research . 45 

3.4  Development  Model  Hardware . 46 

3.5  Dynamic  Analysis . 50 

3.6  Gimbal  Control . 54 

3.7  Difficulties  Encountered  and  other  Considerations . 62 

3.8  Closure . 63 


CHAPTER4 

THE  SIMULATION:  DISCOS . 65 

4.1  Introduction . 65 

4.2  What  is  DISCOS? . 66 

4.3  Dynamic  Analysis . 68 

4.4  Simulation  Results . 69 


5 


CHAPTERS 

EXPERIMENTAL  RESULTS . 76 

5.1  Introduction . 75 

5.2  Defining  Input  Siaping  Effectiveness . 75 

5.3  Open  hoop  Resvdts  and  Discussion . 77 

5.4  Closed  Loop  Frequency  Response . 78 

5.4.1  Experimental  setup . 78 

5.4.2  Shaped  response . 82 

5.4.3  Experimental  insensitivity  curves . 84 

5.4.4  Why  does  it  work  so  well? . 88 

5.4.5  Why  shouldn’t  it  work  so  well? . 89 

5.5  Closed  Loop  Time  Histories . 90 

5.5.1  Slew  profile  and  limitations . 90 

5.5.2  Unshaped  response . 90 

5.5.3  Shap(^  Response . 93 

5.5.3.1  One  mode  shaped . 93 

5.5.3.2  Four  modes  shaped . 96 

5.6  Using  More  Robust  Sequences . 97 

5.7  Identifying  Nonlinearities . 102 

5.8  Using  Different  Controllers . 106 

5.8.1  Breaking  with  impulse  shaping  assumptions . 106 

5.8.2  Designing  different  bandwidth  controllers . 109 

5.8.3  Effectiveness  of  shaping  on  different  controllers . Ill 

5.8.4  Condusions . 115 

5.9  Closure . 116 


CHAPTERS 

CONCLUSIONS  AND  FUTURE  WORK . 119 

6.1  Condusions . 119 

6.2  Suggestions  for  MACE  Project. . 122 

6.3  Future  Work . 125 


REFERENCES . 127 


6 


LIST  OF  FIGURES 

2.1  An  example  of  how  a  step  input  is  convolved  with  an  impulse 

sequence  to  generate  a  ^lap^  command . 24 

2.2  Example  of  posicast  control,  or  a  two-impulse  shaper . 25 

2.3  Three  impulse  sequence  for  (OslOHz,  . 29 

2.4  Insensitivity  curve  for  three  impulse  sequence;  (OslOHz,  . 29 

2.5  An  illustration  of  convolving  a  two  mode  sequence . 31 

2.6  Direct  solution  sequence  for  two  modes,  three  impulses  per  mode 

(2, 20Hz) . 34 

2.7  Convolved  sequence  for  two  modes,  three  impulses  per  mode  (2, 

2OH2) . 34 

2.8  Insensitivity  curve  for  direct  solution  two  mode  sequence 

(D1s2Hz,  Ci=Oi  Ci)2=20Hz,  35 

2.9  Insensitivity  curve  for  convolved  two  mode  sequence  ci>is2Hz, 

Ci=0,  c^=20Hz,  C2=0 .  35 

2. 10  Insensitivity  curve  for  three  impulse  sequence  for  two  modes  at 

9Hz  and  llHz . 39 

3.1  The  MACE  Development  Model . 43 

3.2  MACE  Development  Model  Schematic . 46 

3.3  A  two-axis  giinbal  actuator. . 48 

3.4  Outer  gimbal  stage  to  node  four  accelerometer  open  loop 

gain  transfer  function . 52 

3.5  Outer  gimbal  stage  to  node  four  accelerometer  open  loop 

phase  transfer  function . 52 

3.6  An  illustration  of  the  first  three  vertical  bending 

modes  of  the  DM  model . 53 

3.7  A  block  diagram  of  control  algorithm  for  a  PD  controller  with 

augmented  delay  buffer  states . 60 

3.8  Closed  loop  gain  transfer  function  of  command  to  encoder  angle . 61 

3.9  Closed  loop  phase  transfer  function  of  command  to  encoder 

angle . 61 

4.1  A  schematic  of  what  MACE  looks  like  in  DISCOS . 67 

4.2  DISCOS  model  modal  fii^quendes  in  four  different 

configurations . 68 

4.3  Shai^  and  unshaped  torque  inputs  into  attitude  controller.  1 

mode  (1.65Hz) . 70 

4.4  Bus  inertial  angle  response  to  unshaped  torque  input  into 

attitude  controller. . 70 

4.5  Bus  inertial  angle  response  to  shaped  torque  input  into  attitude 

controller.  1  m^e  (1.65Hz) . 71 


7 


4.6  Payload  angle  response  to  40**  shaped  and  unshaped  slews.  3 

modes  (1.85,  7.5, 15.5Hz) . 72 

4.7  Sensor  deflection  in  y  direction  for  unshaped  40**  slew . 73 

4.8  Sensor  deflection  in  y  direction  for  shaped40*’ dew.  3  modes 

(1.85,  7.5, 15.5H2) . 73 

5.1  Diagram  of  Experimental  Setup . 79 

5.2  Closed  loop  gain  transfer  function  from  commanded  gimbal 

angle  to  vertical  acceleration  at  node  4 . 81 

5.3  Closed  loop  phase  transfer  function  from  commanded  gimbal 

angle  to  vertical  acceleration  at  node  4 . 81 

5.4  Closed  loop  transfer  function  with  1  mode,  3  impulse  shaper 

(1.75Hz) . 83 

5.5  Closed  loop  transfer  function  with  4  mode  direct  solution 

sequence  (1.75, 5.23,  8.68, 9.13Hz) . 83 

5.6  Experimental  insensitivity  curve  for  three  impulse  shaper  for 

first  mode  (1.75H2) . 85 

5.7  Experimental  insensitivity  curve  for  4  mode  direct  solution 

sequence  (1.75, 5.23,  8.68, 9.13Hz) . 85 

5.8  A  comparison  of  experimental  versus  theoretical  insensitivity 

curves  for  one  mode  shaper . 87 

5.9  A  comparison  of  experimental  versus  theoretical  insensitivity 

curves  for  4  mode  ^rect  solution  sequence . 87 

5.10  Profile  of  gimbal  angle  command. . 91 

5.11  Gimbal  ai^e  command  and  response  for  a  smooth  20**  slew. . 91 

5.12  Gimbal  angular  rate  response  to  smooth  20**  slew . 92 

5.13  Node  4  vertical  acceleration  response  to  smooth  20**  slew . 92 

5. 14  Gimbal  angle  command  and  response  for  a  20**  slew  shaped  for 

one  mode  (1.75Hz) . 94 

5.15  Node  4  vertical  acceleration  response  to  20**  slew  shaped  for  one 

mode  (1.75Hz) . 94 

5.16  Fast  Fourier  Transform  of  residual  vibration  at  node  4 

accelerometer  for  20**  slew  shaped  for  1  mode  (1.75Hz) . 95 

5.17  Gimbal  angle  command  and  response  to  20**  slew  shaped  for  4 

modes  (1.75,  5.23, 8.68,  9.13) . 97 

5.18  Node  4  vertical  acceleration  response  to  20**  slew  shaped  for  4 

modes . . . 98 

5. 19  Fast  Fourier  Transform  of  residual  vibration  at  node  4 

accelerometer  for  20**  slew  shaped  for  4  modes . 98 

5.20  A  dose-up  view  of  FFT  of  residual  vibration  amplitude  for  3 

different  impulse  shapers . 101 

5.21  A  dose-up  view  of  FFT  of  residual  vibration  amplitude  for  3 

different  impulse  shapers. . . 101 

5.22  A  comparison  of  FFTs  of  residual  vibration  amplitude  for  four 

different  size  slews,  all  shaped  for  2  modes  (8.68, 9.13Hz) . 104 

5.23  Fast  Foiuier  Transform  of  residual  vibration  amplitude  for  1** 

slew  shaped  for  2  modes  (8.68, 9.13Hz) . 104 

5.24  A  system  diagram  of  MACE . 106 


8 


5.25  Closed  loop  transfer  functions  of  commanded  angle  to  encoder 

output  for  three  different  controllers . Ill 

5.26  Closeup  view  of  FFT  of  residual  vibration  for  lOHz  bandwidth 

controller.  20**  slew  is  shaped  for  2  modes  (8.68, 9.13Hz) . 112 

5.27  Closeup  view  of  FFT  of  residual  vibration  for  3Hz  bandwid^ 

controUer.  20**  slew  is  shaped  for  2  modes  (8.68, 9.13Hz) . 112 

5.28  Closeup  view  of  FFT  of  residual  vibration  for  20Hz  bandwidth 

controller.  20**  slew  is  shaped  for  2  modes  (8.68,  9.13Hz) . 113 

5.29  A  comparison  of  residual  vibration  amplitudes  in  the  6  to  lOHz 

region  for  three  different  controllers . 113 

6.1  A  five  mode  direct  solution  sequence  for  the  MACE  Development 

Model.  (1,75, 4.7, 5.23, 8.68,  9.13Hz) . 124 


9 


LIST  OF  TABLES 


2.1  A  comparison  of  multiple  mode  shaping  methods . 33 

3.1  A  summary  of  the  modal  frequencies  and  damping  ratios  of  the 

MACE  test  artide  in  the  vertical  plane . 47 

3.2  Parameters  used  for  gimbal  modd . 51 

3.3  A  summary  of  the  parameters  for  15Hz  bandwidth  controller 

with  delay  buffer  states . 58 

4.1  DISCOS  model  modal  firequendes  in  four  different 

configurations . 64 

5.1  A  siumnary  of  three  different  bandwidth  controllers . 106 

6.1  Listing  of  impulse  times  and  amplitudes  for  impulse  sequence  in 

figure  6.1 . 119 


11 


12 


Cfiapterl:  Introductim 


18 


INTRODUCTION 

Chapter  1 


1.1  Background 

It’s  a  conunon  problem  that  we  see  all  the  time.  From  bad  suspensions 
that  make  us  carsick  to  the  annoying  blinds  that  flutter  when  the  windows 
are  open,  vibration  is  present  everywhere  in  our  world.  If  we  had  the  ability, 
we  would  all  know  how  to  control  the  vibration  in  these  systems.  Humans 
are  amazingly  good  at  manipulating  things  with  an  uncanny  knowledge  of 
the  vibrational  characteristics  of  structures.  Watching  my  baby  niece,  I 
realized  that  even  she  knew  how  to  shake  her  rattle  to  cause  maximum 
vibration  of  the  beads  inside.  If  she  could  tmderstand  me,  she  could  probably 
shake  the  rattle  such  that  there  would  be  no  vibration  inside  €dso.  The  point 
is  that  humans  have  a  basic  knowledge  of  how  to  make  things  vibrate  and 
how  to  stop  them  from  vibrating.  The  only  remaining  hurdle  is  to  build  the 
tools  to  implement  this  knowledge. 
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With  the  advent  of  robotic  manipulators  in  industrial  applications  and  in 
space,  control  of  vibration  has  become  a  focal  point  of  much  research  in  the 
engineering  community.  For  a  robot  designed  to  do  assembly  in  a 
manufacturing  line,  control  of  vibration  is  critical  in  both  precision  and  in  the 
amount  of  time  required  to  do  its  job.  A  solution  to  the  problem  could  save  a 
company  large  sums  of  money  on  an  annual  basis.  An  example  of  vibration 
problems  in  space  is  the  shuttle  Remote  Manipulator  S3rstem  which  has  its 
fundamental  bending  frequency  in  the  0.5Hz  region.  Having  to  wait  for  this 
vibration  to  settle  after  every  move  is  a  tremendous  loss  of  time  and  money 
for  the  space  program.  As  the  links  in  manipulators  become  lighter  and  more 
flexible  due  to  cost  saving  reasons,  the  need  for  efficient  vibration  control 
schemes  looms  larger. 

In  space  structures,  this  is  becoming  more  and  more  the  case.  With  the 
high  cost  per  unit  of  weight  of  sending  payloads  into  outer  space,  the 
materials  used  in  the  space  structures  tend  to  be  lighter,  and  thus  more 
flexible.  Moreover  there  are  currently  plans  by  NASA  for  building  large 
space  structures.  The  Space  Station  and  the  Earth  Observing  Satellite  are 
two  examples.  The  vibrations  in  these  systems  are  a  problem  on  a  larger 
scale.  The  structural  integrity  becomes  an  issue  when  vibration  is  present. 
Nurre  et  al.  [15]  have  given  a  comprehensive  overview  of  the  issues  involved 
in  controlling  large  spacecraft.  They  start  with  a  definition  of  a  large  flexible 
space  structure  and  go  on  to  survey  some  of  the  modeling,  frequency 
identification,  and  control  issues  in  large  structure  dynamics.  'N^bration 
control  of  large  space  structures  is  the  goal  of  the  work  presented  here. 
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1.2  Input  Shaping 

There  are  two  branches  in  the  field  of  vibration  control  research.  The  first 
is  feedback  control  where  the  algorithm  uses  the  state  of  the  current  S3r8tem 
to  calculate  the  input.  This  is  a  well-researched  field  where  new  ideas  are 
seen  on  a  regular  basis.  The  other  type  of  vibration  control  is  what  I  would 
call  input  shaping.  These  are  generally  open  loop  methods  of  vibratior 
control  that  carefully  manipulate  the  inputs  into  a  system  to  reduce  the 
resulting  vibration.  We  are  interested  in  the  latter  in  this  paper. 


Input  shaping  can  be  further  divided  into  two  subcategories;  those  that 
work  in  real  time,  and  those  that  do  not.  The  latter  type  tends  to  be  more 
computationally  complex  and  have  smoother  profiles.  Many  examples  of  this 
type  of  shaping  are  given  in  the  literature  review  below.  These  methods 
generate  predetermined  inputs  which  are  designed  to  move  a  system  from 
one  point  to  another  without  exciting  vibration.  The  ones  that  work  in  real 
time  tend  to  be  less  computationally  complex  and  can  work  for  any  arbitrary 
input.  Some  examples  of  this  include  lowpass  filters,  digital  notch  filters,  and 
the  subject  of  this  thesis,  impulse  shaping. 

The  most  common,  and  some  may  even  say  the  best,  method  of  input 
shaping  is  to  use  slow  commands.  Usually,  this  will  erase  most  of  the 
vibration  incurred  during  faster  moves.  This  would  be  a  valid  method  alwa3rs 
if  time  was  not  an  issue.  But  it  is.  All  the  attempts  at  input  shaping 
described  below  try  to  minimize  the  time  it  takes  to  move  a  S3rstem  to  a 
certain  point  without  leaving  any  reddual  vibration.  An  input  shaping 
algorithm  always  has  to  be  measured  against  the  cost  in  time  delay  as  well  as 
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its  efifectiveness.  The  bottom  line  for  an  input  shaper  is  the  time  it  takes  for 
the  system  to  settle  down  to  an  acceptable  vibration  levd. 


1.3  Literature  Review 

Some  of  the  early  work  in  input  shaping  includes  Swigert  [23],  who 
showed  ways  to  calculate  torques  on  simple  rotational  structures  which 
resulted  in  zero  residual  vibration.  He  started  with  a  simple  state  space 
model  and  found  waveforms  that  met  the  zero  vibration  terminal  boundary 
conditions.  He  also  included  a  penalty  on  residual  amplitude  that  allowed 
limited  errors  in  plant  frequency  identification.  Aspinwall  [2]  used  sine 
series  to  create  forcing  functions  that  caned  vibration  in  single  degree  of 
fii^edom  undamped  systems.  By  using  a  finite  Fourier  series  expansion,  he 
was  able  to  find  the  coefficients  to  the  sine  series  that  minimize  the  fiwquenqy 
content  over  a  wide  range  of  firequendes.  This  was  an  effective  method  but 
was  only  considered  for  undamped  s^^tems  and  the  forcing  functions  were 
lengthy.  Farrenkopf  [8]  used  Calculus  of  Variations  to  solve  the  optimal 
slewing  profile  problem  with  initial  and  final  boundary  conditions  set.  He 
found  that  on  a  simple  one  mode  system,  it  is  difficult  to  not  exdte  the 
structural  mode  when  at  <  5,  where  a  is  the  modal  frequeni^  and  t  is  the 
total  time  of  Ike  move. 

Meckl  [12]  did  some  similar  work  in  generating  smooth  force  profiles  that 
are  shaped  to  reduce  energy  levels  at  the  modal  firequendes.  He  used  the 
sum  of  ramped  sinusoids  and  chose  the  coeffidents  by  minimizing  the 
spectral  magnitudes  over  a  range  of  firequendes  and  the  time  of  the  move.  Bte 
also  developed  profiles  that  are  fundamental  versine  curves  which  can  be 
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applied  to  velocity  limited  systems.  A  closed  loop  implementation  of  these 
methods  was  described  where  the  force  profiles  are  double  integrated  to 
generate  reference  position  commands  and  fed  forward. 

Eisler,  Segalman»  and  Robinett  [6]  have  used  recursive  quadratic 
programming  and  a  finite  element  model  to  generate  approximate  minimum 
time  trajectories  for  a  two  link  flexible  robot.  Boundary  conditions  were  set 
on  the  mid  position,  energy,  straight-line  tracking,  and  motor  torque  limits  to 
optimize  the  parameters  of  the  move.  Their  results  showed  good  accuracy  in 
tip  tracking,  and  it  was  found  interestin^y  that  the  calculated  torque  profiles 
were  similar  even  for  large  variances  in  plant  configuration  and  modal 
fi:«quencie8.  Turner  and  Junkins  [25]  have  looked  at  optimal  single-axis 
reorientations  of  a  structure  with  flexible  arms.  Necessary  conditions 
stemming  from  Pontryagin's  principle  were  applied  to  arrive  at  the  move 
profiles.  Also,  an  extension  of  the  algorithm  to  include  kinematic 
nonlinearities  was  presented. 

Inverse  dynamics  is  another  method  of  input  shaping  that  has  been 
researched.  Many  recent  textbooks  describe  the  method  of  inverting  the 
system  model  to  calculate  the  required  inputs  for  a  given  desired  trajectory. 
Asada,  Ma,  and  Tokumaru  [1]  present  a  technique  that  reduces  the 
computation  complexity.  They  use  what  they  term  a  virtual  rigid  Knk 
coordinate  system  to  simplify  the  boundary  conditions  that  need  to  be  set. 

Prucz,  Soong,  and  Reinhom  [18]  present  an  interesting  twist  on  the  use  of 
impulses  in  vibration  control.  In  ^Pulse  Control,”  inq>ul8es  are  input  into  the 
S3r8tem  to  limit  the  system  response  to  a  certain  threshold  level  at  all  times. 
They  propose  an  efficient  algorithm  to  calculate  the  optimal  times  and 
amplitudes  of  these  impulses  to  keep  the  vibration  level  low. 
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The  pioneer  of  the  impulse  shaping  method  may  arguably  have  been 
Smith  [22].  Way  back  in  195S,  he  proposed  what  he  called  *Posicast  Control.” 
It  requires  that  a  step  input  be  divided  into  two  equal  parts,  the  second  half 
being  delayed  in  time  by  half  the  period  of  the  modal  frequency.  By  doing 
this,  the  plant  can  be  moved  without  ezdting  its  resonant  frequency.  He  also 
proposed  a  closed  loop  implementation  of  posicast  control  that  could  work  to 
reject  distuihances  in  the  system.  This,  however,  reqiures  an  exact  model  of 
the  system  and  is  not  very  robust  to  errors  in  system  identification. 

Singer  [20]  wrote  the  first  authoritative  work  on  the  modem  technology 
referred  to  as  impulse  shaping.  He  discovered  a  way  to  use  the  mathematics 
behind  the  response  of  second  order  systems  to  impulses  to  create  real  time 
filters  that  produce  vibration  reducing  inputs.  By  setting  the  residual 
vibration  amplitude  to  go  to  zero,  he  was  able  to  generate  sequences  of 
impulses  that  did  not  excite  the  modal  frequencies.  Real  time 
implementation  was  achieved  through  convolution  of  these  impulses  with  any 
arbitrary  input.  The  breakthrough  in  this  technology  was  the  ability  to 
operate  in  the  presence  of  errors  in  system  parameter  identification.  He 
derived  his  method  in  both  the  firequency  and  time  domains  and  presented 
experimental  results  that  showed  the  method  to  be  extremely  effective  in 
simple  structures.  He  also  did  a  thorough  comparison  study  of  the  method  to 
digital  notch  filters.  He  found  that  the  impulse  shaping  method  was  superior 
in  decreasing  the  amount  of  vibration  and  was  also  more  robust  to  system 
parameter  variations. 

Hyde  [10]  introduced  a  method  of  calculating  multiple  mode  impulse 
sequences  throui^  a  direct  solution  method  which  decreases  the  total  time 
length  of  Singer’s  sequences.  Though  more  complex  in  computation,  these 
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sequences  had  fewer  impulses  than  the  previous  multiple  mode  sequences 
and  had  a  time  savings  of  20-30%.  He  verified  his  restilts  in  both  linear  and 
non-linear  simulations.  Singhose  [21]  also  proposed  an  extension  to  Singer’s 
original  method.  He  showed  that  more  robustness  can  be  obtained  by 
specifying  a  tolerance  level  in  vibration  amplitude.  By  releasing  the  zero 
vibration  (x>nstraint,  a  larger  region  of  insensitivity  can  be  achieved  by 
modifying  the  basic  three  impulse  sequence  for  one  mode  developed  by 
Singer. 

There  have  been  a  handful  of  papers  written  since  the  advent  of  impulse 
shaping  by  people  who  have  successfully  implemented  the  method.  Christian 
[4]  built  a  two  link  flexible  robot.,  called  the  ‘Tlezbot,’*  which  was  specifically 
designed  to  study  vibration  control.  He  compared  impulse  shaping  to  slow 
acceleration  moves  and  showed  that  impulse  filtering  is  a  more  effective  way 
to  reduce  vibration  in  systems  with  constant  vibrational  frequencies.  Tzes, 
Englehart,  and  Yurkovich  [26,  27]  have  experimented  with  an  adaptive  form 
of  impulse  shaping  that  uses  firequency  domain  information  firom  the  plant  to 
control  systems  with  time  varying  frequencies.  They  use  a  real  time 
frequency  identification  scheme  to  compute  modal  firequendes  and  update  the 
spacing  between  impulses.  By  doing  this,  they  were  able  to  use  the  shorter 
two  impulse  sequences  and  achieved  good  robustness  to  variances  in  plant 
frequences.  Murphy  and  Watanabe  [15]  offer  a  digital  analysis  of  the 
impulse  shaping  formulation  and  show  that  the  digitized  form  is  a  class  of 
notch  filters  that  place  zeroes  to  cancel  the  plant  modal  roots.  They  also 
propose  an  arbitrary  sampling  rate  digital  shaping  filter  that  allows  the 
implementation  of  these  filters  in  varying  frequency  systems. 
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1.4  Overview  of  Thesis 

The  goal  of  my  work  is  to  answer  some  of  the  implementation  issues 
involved  in  impulse  shaping.  Most  of  the  past  researdi  has  concentrated  on 
theory  and  application  to  simple  systems.  Simulations  have  been  done  on 
more  complex  systems,  but  there  still  isn’t  enough  experimental  data  to  raise 
the  necessary  questions  of  implementation.  A  theoretician  can  be  very 
confident  about  her  design  on  paper,  but  there  will  always  be  unmodeled 
factors  that  arise  when  working  in  real  hardware.  This  thesis  will  answer 
some  of  the  questions  that  have  to  do  with  implementing  impulse  shaping  on 
a  complex  structure. 

A  comprehensive  overview  of  the  impulse  shaping  theory  is  given  in 
Chapter  2.  I  will  go  over  Singer’s  equations  and  then  cover  the  modifications 
made  by  Hyde  and  Singhose.  Then,  I  will  present  a  comparison  of  the 
different  methods  of  computing  multiple  mode  sequences.  Some  tradeofib  will 
be  shown  between  Singer's  convolution  method  and  Hyde’s  direct  solution 
method.  I  will  end  the  chapter  with  some  considerations  in  choosing  the  most 
efficient  sequence  for  a  multiple  mode  system. 

Chapter  3  will  reveal  the  testbed  on  which  the  experiments  will  take 
place.  The  Middeck  Active  Control  Experiment  (MACE)  is  a  planned  space 
shuttle  mission  project  that  will  be  shown  to  be  an  ideal  platform  on  which  to 
generate  valuable  data  about  the  effectiveness  of  impulse  shaping  on 
complicated  structures.  I  will  describe  the  current  hardware  for  the  project 
and  detail  some  control  algorithms  for  the  gimbal  actuator. 


Nonlinear  simulation  will  be  the  topic  of  Chapter  4.  DISCOS  (Dynamics 
Interaction  Simulation  of  Controls  and  Structures)  will  be  used  as  the  tool 
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for  modeling  and  simulating  an  early  version  of  the  MACE  hardware. 
Impulse  shaping  results  for  both  open  loop  and  closed  loop  runs  will  be 
shown. 

Experimental  results  from  the  hardware  will  be  presented  in  Chapter  5. 
A  description  of  the  flight  objectives  for  impulse  shaping  will  be  explained  to 
define  the  problem.  In  both  frequency  domain  results  and  time  domain 
results,  impulse  shaping  will  be  shown  to  be  very  effective  for  amall  ang^e 
moves.  Insensitivity  curves  will  be  compared  from  experimental  to 
theoretical  data.  A  further  analysis  will  be  presented  concerning  the  non- 
linearities  associated  with  large  angle  movements.  Also,  the  involvement  of 
the  controller  bandwidth  in  the  effectiv^n^  of  shapers  will  be  shown.  I  will 
show  through  these  experiments  that  although  impulse  shaping  can  work 
very  well  in  certain  configurations,  it  does  have  its  shortcomings  as  non- 
linearities  arise. 

Finally,  in  Chapter  6,  I  will  reach  some  conclusions  about  all  the  data 
presented  emd  look  to  the  future  with  suggestions  for  both  the  MACE  project 
and  further  research  on  impulse  sha|nx^. 
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THE  TECHNOLOGY: 
IMPULSE  SHAPING 


Chapter  2 


2«1  Introduction 

This  chapter  deals  ^th  the  theoretical  basis  of  impulse  shaping.  The 
objective  of  this  type  of  input  shaping  is  to  use  a  series  of  impulses  to 
suppress  the  residual  vibration  of  a  dynamic  S3rstem.  By  solving  for  a  set  of 
vibration  constraint  equations  based  on  the  response  of  systems  to  impulses, 
a  sequence  of  unique  impiilses  can  be  formulated  which,  when  convolved 
with  the  input  to  a  system,  will  suppress  the  resonant  mode  vibration.  The 
following  sections  are  meant  to  give  a  comprehensive  overview  of  the 
algorithm  for  generating  the  sequences  and  to  lend  some  insight  into  the 
method.  The  authors  cited  should  be  consulted  for  more  detailed 
explanations. 


2.2  Implementing  the  Sequences 


One  first  needs  to  understand  how  these  impulse  sequences  are  used. 
Impulses  are  not  actually  input  into  the  system,  as  that  would  be  hard  on 
the  amplifiers,  and  it  would  not  achieve  any  desired  move.  Instead,  the 
impulse  sequence  is  used  as  a  filter  to  “shape”  an  arbitrary  move.  This  is 
accomplished  by  performing  a  convolution  of  the  original  input  with  the 
impulse  sequence.  By  doing  this,  we  are  projecting  the  vibration  reducing 
attributes  of  the  impulse  sequence  onto  an  arbitrary  move.  This  presents 
one  of  the  many  advantages  of  imptilse  shaping,  which  is  that  it  can  be  used 
in  real  time,  with  any  type  of  input.  An  example  of  what  a  “shaped”  input 
looks  like  is  shown  in  figure  2.1. 


rnmmanrf  Thrcc  impulsc  filter  Shaped  Command 
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FIGURE  2.1:  An  example  of  how  a  step  input  is  convolved  with  an  impulse  sequence  to 

generase  a  shaped  command. 

2.3  Calculating  Single  Mode  Sequences 

Much  like  Smith's  posicast  control,  impulse  shaping  delays  parts  of  the 
original  input  based  on  the  natural  frequency  of  the  plant.  A  good  way  to 
visualize  what  is  happening  is  shown  in  figure  2.2.  A  simple  pendulum  on 


a  string  is  chosen  as  the  plant  to  be  controlled.  The  desired  move  is  to  go 
firom  point  A  to  point  B  without  causing  the  mass  to  swing  at  the  end  of  the 
move.  The  pendulum  motion  is  the  resonant  vibration  of  this  system.  The 
move  is  split  into  two  parts.  The  pendulum  is  first  moved  half  the  distance 
to  point  B.  Then,  after  the  mass  has  swung  all  the  way  to  the  other  side,  the 
other  half  of  the  move  is  completed.  The  desired  move  has  been  achieved 
without  any  residual  vibration.  It  is  a  very  simple  idea  but  it  is  the 
foundation  for  a  very  sophisticated  vibration  control  strategy  that  has  many 
important  applications. 


FIGURE  2.2:  Example  of  posicast  ccmtrol,  or  a  two-impulse  shaper. 


The  modem  technology  that  we  call  impulse  shaping  is  a  few 
generations  removed  firom  Smith's  posicast  control  idea.  The  original  work 
done  in  this  area  by  Singer  starts  with  the  assumption  that  most  dynamic 
systems  can  be  modeled  as  a  stun  of  second  order  linear  systems.  The 
response  of  a  second  order  system  to  an  impulse  is  governed  by  the 
following  equation: 


y(t)  =  A  sin  (©  V 1-C^  (t-t©)) , 


(2.1) 
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where  o  is  the  tindamped  natural  frequency  of  the  system,  ^  is  the  damping 
ratio  of  the  system,  i^is  the  time  of  the  impulse,  and  A  is  the  amplitude  of 
the  impulse.  By  superposition,  we  can  then  conclude  that  the  response  to  a 
series  of  impulses  will  be  the  sum  of  the  responses  to  each  of  the  impulses. 
We  can  then  write: 


y(t)  =  X  sin  (co  V 1-C^  (t-to) , 


(2.2) 


i=l 


where  n  is  the  number  of  impulses.  From  this  we  can  derive  an  equation 
for  the  amplitude  of  vibration  at  the  time  of  the  last  impulse.  This  equation 
describes  the  residual  vibration  amplitude  envelope. 


Amp 


.V( 


X  Bi  cos  +  (X 

i*l  /  \i»l 

Bi  =  AiC-C  «>(«--“) 

(|)i  =  (0  V  1  -  ti 


(2.3) 


To  obtain  a  vibration-firee  impulse  sequence,  we  want  this  expression  to 
equal  zero.  This  requires  that  both  squared  terms  in  equation  2.3 
independently  go  to  zero,  yielding  the  following  constraint  equations: 

^  Ai  sin  (©  (t-to)  =  0 

X  Ai  cos  (o)  V 1-C^  (t-ti))  =  0 


(2.4) 
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Given  that  we  know  o)  and  the  undetermined  variables  are  the 
impulse  amplitudes  Ai  and  the  impulse  times  ti.  More  constraint  equations 
are  added.  We  realize  that  only  the  relative  amplitudes  are  important  and 
that  the  starting  time  of  the  impulse  sequence  is  arbitrary,  so  we  list  two 
more  constraints: 


(2.5) 


The  time  of  the  first  impulse  is  set  to  be  zero  and  the  amplitudes  are 
normalized  to  maintain  the  same  final  stc  iy  state  value  of  the  original 
input  during  implementation.  We  now  hav^  four  equations  and,  since  each 
impulse  has  two  variables  (time  and  amplitude),  two  impulses  can  be  found 
that  will  satisfy  these  equations.  The  resulting  two  impulse  sequence,  when 
implemented  with  a  step  input,  is  not  unlike  Smith's  posicast  control. 

The  breakthrough  in  impulse  shaping  is  its  robustness  factor.  By 
robustness  we  mean  that  the  method  should  be  fairly  effective  even  for 
slight  errors  in  system  parameter  identification  (o  and  C)-  This  robustness 
constraint  is  modeled  by  taking  the  first  derivative  of  equation  2.4  with 
respect  to  oi  and  setting  it  to  zero.  Singer  proves  that  the  derivative  of  2.4 
with  respect  to  o>  is  the  same  as  the  derivative  with  respect  to  C-  So  the 
following  derivative  constraints  require  that  the  impulse  train  have  little 
residual  vibration  for  errors  in  both  the  system  natural  frequency  and 
damping  ratio. 


(2.6) 


^  Ai  ti  sin  (ti  00  V 1-C^)  =  0 

i-l 

AitiC-W^  cos  (tiwV  1-C^)  =  0 

i>l 


We  now  have  six  equations  which  can  be  satisfied  by  three  impulses. 
This  is  a  basic  three  impulse  sequence  for  one  mode  that  is  fairly  robust  to 
system  parameter  uncertainty.  To  gain  more  robustness,  we  can  add 
another  constraint  by  taking  a  second  derivative,  yielding  a  fourth  impulse. 
We  can  actually  go  on  and  on  by  adding  more  and  more  derivative 
constraints.  The  tradeoff,  however,  is  the  length  of  the  impulse  train.  Every 
additional  derivative  constraint  adds  an  additional  length  to  the  sequmice 
that  is  equal  to  half  the  period  of  the  firequency  being  controlled.  The  longer 
the  impulse  train,  the  longer  the  delay  will  be  in  implementation.  An 
example  of  a  three  impulse  sequence  for  a  lOHz  resonant  system  is  shown 
in  figure  2.3.  Figure  2.4  shows  an  insensitivity  curve  for  the  three  impulse 
sequence.  The  insensitivity  curve  is  way  to  graphically  portray  the 
robustness  of  an  impulse  sequence  by  showing  the  residual  vibration 
amplitude  as  a  function  of  frequen^r.  The  curve  is  a  calculation  of  the 
residual  vibration  amplitude  response  for  each  frequency  to  the  given 
impulse  sequence,  and  plotted  on  a  linear  scale  as  a  percentage  of  the 
residual  vibration  amplitude  for  a  unit  impulse.  The  plot  is  thus  a 
measurement  of  how  much  vibration  amplitude  a  certain  impulse 
sequence  can  cancel  at  each  firequency  for  a  given  range. 


HMECSEQ 


Three  impulse  sequrace  for  ooslOHz, 
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Insensitivity  curve  for  diree  inqmlse  sequence;  ooslOHz, 
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2.4  Calculating  Multiple  Mode  Sequences 

2.4.1  Convolution 

This  algorithm  is  now  ready  to  be  extended  to  multiple  mode  problems. 
The  method  suggested  by  Singer  is  to  convolve  the  impulse  sequences  for 
each  of  the  individual  modes.  Mathematically,  this  involves  performing  an 
outer  sum  of  the  impulse  times  and  an  outer  product  of  the  impulse 

m 

amplitudes.  The  resulting  sequence  has  Jl  impulses,  and  a  length  of 

m 

X  ht  where  m  is  the  ntunber  of  modes,  n  is  the  number  of  impulses  for  each 

mode,  and  /  is  the  length  of  each  sequence.  An  easier  way  to  think  of 
convolution  is  to  picture  a  one  mode  impulse  shaper  as  a  black  box  filter 
which  splits  an  input  into  three  parts,  delaying  the  latter  two.  Convolution 
puts  these  black  boxes  in  series  such  that  the  "filtered”  output  from  the  first 
shaper  is  used  as  input  into  the  second  shaper  and  so  forth.  The  resulting 
impulse  sequence  is  the  output  from  the  last  shaper  in  the  series.  Figure 

2.5  shows  the  type  of  sequence  that  would  result  firom  convolving  two  3- 
impulse  shapers. 

Convolved  sequences  are  easy  to  calculate  since  they  involve  only  finding 
the  sequences  for  one  mode  problems  and  some  simple  math  in  joining 
them.  They  are  also  easy  to  implement  in  computer  code  by  using  this  black 
box  idea  and  programming  the  individual  one  mode  shapers  in  series.  This 
results  in  less  code  and  also  less  calculations  for  the  computer,  whidi  is  an 
advantage  when  systems  are  working  at  high  servo  rates.  The 
disadvantage  to  convolving  is  the  length  of  the  resulting  sequence.  As  will 
be  shown  in  the  next  section,  shorter  sequences  can  be  found  to  satisfy  the 
same  vibration  constraint  equations. 
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FIGURE  2.5;  An  iUustratioa  of  ctxivolving  a  two  mode  sequence 


2.4^  Direct  Solution 

Hyde  has  detailed  a  direct  solution  method  for  finding  impulse 
sequences  for  multiple  mode  problems.  The  algorithm  basically  lays  out  the 
vibration  constraint  equations  for  all  the  modes  and  solves  them 
simultaneously.  For  example,  if  we  are  solving  for  3  impulse  sequences,  we 
can  generate  the  constraint  equations  for  an  arbitrary  number  of  modes  m; 


^  Ai  sin  (ofti  Vl-(^  (tn-to)  =  0 

i-l 

^  Ai  cos  (q^  VT^  (t„-ti)) »  0 

1-1 

^  Ai  ti  sin  (osj  V  1-Cj  (tn-to)  =  0 

1-1 

Ai  ti  cos  (co^  V  1-Cj  (tn-ti))  =  0 

(2-6) 

ti  *  0 

tA,  =  1 
^1 

where  the  top  four  equations  are  repeated  for  each  mode  J.  We  now  have  a 
system  of  4in  -i*  2  equations,  which  is  satisfied  by  2m*i>  1  impulses.  So  for  a 
two  mode  case,  the  direct  solution  method  results  in  five  impulses  where 
convolution  would  come  up  with  nine.  The  calculated  sequence  not  only 
has  fewer  impulses  than  the  convolved  sequence,  but  also  is  shorter  in 
length.  The  savings  in  length  is  not  always  predictable  but  experience  has 
shown  that,  for  three  impulse  sequence  modes  of  less  than  five  modes,  the 
savings  can  be  up  to  30%.  When  we’re  dealing  with  slow  modes,  the 
savings  in  time  can  be  a  big  factor. 

The  method  does  have  its  shortcomings,  however.  So  fiur,  an  easy  way  to 
solve  the  nonlinear  constraint  equations  has  not  been  found.  Currently, 
linear  approximations  generated  <hrou|^  optimization  are  used  for  non¬ 
linear  solving  routines.  The  procedure  is  complex,  time-consuming,  and 
does  not  always  find  an  exact  solution.  The  author  has  many  times 
experienced  the  frustration  of  spending  much  time  for  a  multiple-mode 
problem  only  to  find  that  the  solver  does  not  converge.  Contributions  firom 
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ciirrent  research  on  developing  closed-form  solutions  to  multiple  mode 
problems  could  make  this  method  much  more  attractive.  The  robustness  of 
the  convolution  method  is  lost  also  in  the  direct  solutions  as  is  shown  in  the 
following  section. 


2.4^  Convolution  vs.  Direct  Solution 

The  two  approaches  to  solving  multiple  mode  problems  produce 
somewhat  different  insensitivity  curves.  Although  they  satisfy  the  same 
zero  vibration  constraint  equations,  the  robustness  of  the  two  types  of 
sequences  varies  considerably.  Insensitivity  curves  are  usually  plotted  as  a 
function  of  a  normalized  frequency.  This  means  the  robustness  is  a 
function  of  a  percentage  of  the  frequency,  not  absolute  frequency  (Hz,  for 
example).  So  a  IHz  mode  sequence  may  have  5%  residual  vibration 
robustness  between  .85Hz  and  1.15Hz,  but  a  lOHz  mode  sequence  will  have 
the  same  robustness  between  8.5Hz  and  11.5Hz,  an  interval  ten  times 
greater  than  that  of  the  IHz  sequence.  So,  in  absolute  terms,  we  can  say 
that  the  lOHz  sequence  is  ten  times  more  robust  than  the  IHz  sequence. 
This  is  a  fortunate  feature  of  impulse  shaping  since  the  higher  modes  are 
generally  harder  to  identify  precisely  and  they  tend  to  move  around  more. 

In  comparing  the  insensitivity  curves  for  convolved  and  direct  solution 
sequences,  we  notice  a  vast  difference  in  robustness  for  the  higher  modes. 
The  convolved  sequence  maintains  the  same  robustness  in  the  multiple 
mode  sequence  as  in  the  one  mode  sequence.  This  means  that  if  we 
compare  the  insensitivity  curve  in  the  area  of  20Hz  for  a  one  mode  sequence 
for  20Hz  and  a  two  mode  convolved  sequence  for  lOHz  and  20Hz,  they  are  the 
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FIGURE  2.6:  Direct  solutkn  sequence  for  two  modes,  three 
inqMilses  per  mode  (2, 20Hz). 


FIGURE  2.7:  Convolved  sequence  for  two  modes,  three 
impulses  per  mode  (2, 20Hz). 


VTORATION  AMPLITUDE 


FIGURE  2.8:  In^nsitivity  curve  for  diiect  solutirai  two  mode  sequence; 
Qi)i=:2Hz  ^2=0. 


FIGURE  2.9:  Insensitivity  curve  for  convolved  two  mode  sequence; 
cdi=2Hz  eo2“20Hz  ^2=0. 


same.  However,  in  the  direct  solution  sequences,  the  higher  modes  seem  to 
inherit  the  robustness  curve  shape  of  the  slowest  mode  that  is  shaped.  So  if 
we  compare  the  same  two  sequences  solved  by  the  direct  solution  method, 
we  find  that  the  insensitivity  curve  in  the  area  of  20Hz  for  the  two  mode  case 
is  only  as  robust  as  the  lOHz  region,  about  half  as  much  as  the  convolved 
sequence.  Figures  2.6  and  2.7  show  the  impulse  sequences  using  both  the 
convolution  and  direct  solution  method  for  a  two  mode  problem  using  three 
impulses  per  mode.  The  two  modes  are  placed  an  order  of  magnitude  apart 
to  magnify  the  difference  between  the  two  methods.  Figures  2.8  and  2.9 
show  the  corresponding  insensitivity  curves  for  the  two  sequences.  Notice 
the  difference  in  insensitivities  for  both  sequences  in  the  20Hz  region. 

A  peculiar  feature  of  this  combination  of  modes  is  that  the  direct 
solution  method  yields  a  longer  sequence  than  the  convolved  sequence.  This 
has  to  do  with  the  fact  that  the  ratio  of  the  two  frequencies  is  rather  large. 
There  is  no  current  method  to  predict  when  this  will  happen,  but  we  do 
know  that  for  a  two  mode  sequence,  a  ratio  of  frequencies  that  is  greater 
than  3  will  yield  a  longer  direct  solution  sequence.  Further  investigation 
into  closed  form  solutions  for  the  direct  solution  method  may  reveal  what 
causes  this  to  happen.  For  most  cases,  however,  the  direct  solution  method 
displays  a  significant  savings  in  length. 

To  summarize,  we  can  say  that  there  are  certain  advantages  and 
disadvantages  to  each  method.  Convolution  is  simple  to  calculate,  easy  to 
program  and  implement,  and  maintains  high  robustness.  Direct  solution 
produces  a  shorter  sequence,  but  is  more  complex  in  calculation,  harder  to 
implement,  and  loses  some  robustness.  Different  applications  will  call  for 
different  methods.  When  reducing  the  time  lag  in  the  move  is  critical. 
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direct  solution  is  the  best  choice.  An  application  with  hi^  modes  which 
tend  to  move  around  a  lot  will  call  for  convolution.  And  in  many  cases,  a 
combination  of  the  two  is  ideal.  A  list  of  the  qualities  of  both  methods  is 
given  in  table  2.1. 


TABLE  2.1:  A  cott^arison  of  multiple  mode  shq)ing  methods, 
m  is  the  number  of  modes,  n  is  die  number  of  impidses  per  mode, 
and  /  is  the  length  of  each  sequence. 


2.5  Finding  the  Most  Efficient  Sequence 

There  needs  to  be  some  sort  of  strategy  in  finding  the  “best”  sequence  for 
a  particular  application.  Ideally  what  we  would  like  to  do  is  to  somehow 
minimize  the  length  of  the  sequence  while  maximizing  the  robustness  of 
the  insensitivity  curve.  There  are  some  variables  to  consider.  We’ve 
already  considered  direct  solution  versus  convolution.  The  other  two 
variables  are  the  number  of  modes  to  shape,  and  the  number  of  derivative 
constraints  to  use 


In  determining  the  number  of  modes  to  be  shaped,  the  first  thing  to 
realize  is  the  repeating  nature  of  impulse  shaping.  Because  we  are  using 
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impxilses  at  half  periods  of  the  frequency,  we  also  happen  to  cancel  modes 
which  are  odd  harmonics  of  the  original  frequency.  By  sheer  coincidence, 
this  can  often  eliminate  the  need  to  shape  for  certain  higher  modes, 
resulting  in  shorter,  less  complicated  sequences. 

Also,  as  mentioned  earlier,  higher  mode  sequences  are  much  more 
robust  than  lower  mode  sequences.  This  can  also  be  used  in  our  favor.  If 
two  fairly  frst  resonant  peaks  are  dose  together,  usually,  a  sequence  for  one 
of  the  modes,  or  a  sequence  for  the  frequency  in  between  the  two  modes  will 
suffice.  The  higher  the  modes,  the  more  this  applies.  These  are  two 
inherent  features  of  impulse  shaping  that  we  can  use  in  our  favor  to 
simplify  our  algorithm. 

The  addition  of  a  derivative  constraint  adds  one  more  impulse  to  the  tail 
end  of  a  one  mode  sequence.  It  also  tends  to  lengthen  the  sequence  for 
multiple  mode  problems.  In  the  recent  research  done  on  impulse  shaping, 
a  three  impulse  sequence  has  been  found  to  be  adequate  for  most 
applications.  However,  the  number  of  derivative  constraints  for  a  mode 
should  depend  on  the  uncertainty  of  the  plant  or  the  tendency  of  that  mode 
to  change. 

We  can  also  use  less  derivative  constraints  when  modes  are  dose  to  one 
another.  Singhose  has  proposed  a  method  which  takes  advantage  of  modes 
which  are  close  together  to  yield  shorter  sequences.  His  method  is  actually 
equivalent  to  a  direct  solution  sequence  for  two  modes  without  derivative 
constraints.  The  argument  centers  around  the  fact  than  when  2  shaped 
modes  are  close  together,  their  robustaiess  tend  to  “overlap.”  Whereas  a 
single  two  impulse  shaper  is  not  very  robust,  if  two  2*impulse  sequences  are 
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placed  close  to  each  other  in  frequency,  the  region  in  between  the  two  modes 
is  still  pretty  robust.  Given  a  certain  tolerance  level,  the  two  modes  can  be 
placed  an  optimal  distance  apart  to  gain  maximum  robustness.  This  is 
used  in  place  of  three>impulse  sequences  for  both  modes  or  one  three- 
impulse  sequence  for  the  frequency  in  between  the  two  modes.  The 
resulting  insensitivity  curve  may  look  something  like  figure  2.7.  It  shows  a 
three  impulse  sequence  that  has  been  solved  for  two  modes  around  lOHz 
(9Hz  and  llHz).  The  figure  may  be  compared  against  a  similar  three 
impulse  sequence  for  the  lOHz  mode  alone. 


FIGURE  2.10:  Insensitivity  curve  for  three  impulse  sequence  for  two 
modes  at  9Hz  and  1  IHz. 


We  have  covered  three  variables  we  can  control  to  come  up  with  the  most 
efficient  impulse  sequence:  convolution  vs.  direct  solution,  the  number  of 
modes  to  shape,  and  the  number  of  derivative  constraints.  These  variables 
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can  be  used  in  our  favor  to  optimize  the  tradeoff  between  robustness  and 
time  delay.  The  robustness  &ctor  and  the  ability  to  implement  the  shapers 
in  real  time  are  the  main  advantages  of  impulse  shaping  technology.  The 
tools  are  now  in  place  to  exploit  these  factors  to  produce  more  efficient 
sequences. 


THE  TESTBED:  MACE 

Chapter  3 


3.1  Introduction 

One  of  the  most  suitable  applications  of  impulse  shaping  is  in  space 
structures.  Residual  vibration  is  a  common  problem  in  large  flexible  space- 
borne  systems  and  has  been  the  focus  of  extensive  recent  research  effort  in 
the  aerospace  and  mechanical  engineering  communities.  For  this  reason, 
the  Middeck  Active  Control  Experiment  (MACE)  has  been  chosen  as  the 
testbed  for  verification  of  the  issues  concerning  impulse  shaping  presented 
in  this  work.  This  chapter  will  describe  the  MACE  project  and  detail  some 
of  the  work  done  in  control  derivation. 

3*2  Project  Motivation  and  Agenda 


NASA  initiated  efforts  in  1988  to  investigate  the  behavior  of  on-orbit 
closed  loop  dynamic  systems  which  utilize  Controlled  Structures 
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Technology  (CST).  CST  refers  to  the  use  of  active  control  techniques  in 
highly  flexible  structures  with  closely  coupled  modes  within  the  controller 
bandwidth.  One  of  the  products  of  this  effort  was  the  MACE  project,  a 
proposed  Space  Shuttle  middeck  experiment.  The  goal  of  the  project  was  to 
provide  a  platform  for  analyzing  closed  loop  behavior  in  a  zero  gravity 
environment.  This  project  is  on-going  at  the  MIT  Space  Engineering 
Research  Center  and  managed  by  the  NASA  Langley  Research  Center. 
The  latmch  date  is  currently  set  for  the  summer  of  1994. 

A  multibody  platform  was  chosen  as  the  template  for  the  MACE  design. 
It  was  desired  that  the  results  of  the  experiment  be  useful  to  planned  future 
space  missions.  The  design,  therefore,  takes  after  a  scaled  down  version  of 
a  larger  functional  proposed  design  such  as  the  Earth  Observing  System 
(EOS),  the  Geostationary  Platform  (GEOS),  or  other  exploratory  or  military 
designs.  The  multibody  configuration  is  conducive  to  studying  many  issues 
associated  with  active  control: 

*  Three  dimensional  behavior 

*  Disturbance  isolation 

*  Time  varying  d3niamics 

*  Multiple  interacting  control  systems 

These  featrires  allow  MACE  to  be  the  subject  of  various  types  of  analysis 
whose  findings  can  be  useful  in  futiire  space  flights. 

The  experiments  will  be  run  on  the  Space  Shuttle  middeck  because  of  the 
low  relative  cost  of  running  experiments  there  and  the  ease  of  access  to  the 
facility.  The  middeck  dimensions  are  mission  specific  because  of 
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miscellaneous  hardware  being  stored  in  the  area,  but  it  can  most  certainly 
accommodate  a  test  article  76*  (1.9m)  long.  Maximum  height  of  the 
middeck  is  about  88*  and  maximuni  depth  is  about  60”.  This  restricted  the 
outer  dimensions  of  MACE.  The  flight  plan  also  calls  for  the  MACE 
hardware  being  stowed  in  a  midde<^  locker  with  about  two  cubic  feet  of 
volume.  This  necessitated  the  need  for  a  modular  design  such  that  the 
structure  could  be  taken  apart  and  reassembled  easily  by  the  shuttle  crew. 
Other  restrictions  on  power,  weight,  and  safety  dictated  much  of  the  design. 
The  MACE  hardware  is  shown  in  figure  3.1. 

The  main  science  objective  is  to  investigate  the  extent  to  which  closed 
loop  behavior  of  an  on-orbit  dynamic  system  can  be  predicted  on  earth.  It  is 
impossible  to  reasonably  simulate  on  the  ground  the  environment  that  the 
test  article  would  see  in  space.  The  prediction,  thus,  must  come  from  a 
combination  of  theoretical  modeling  and  test  results  in  a  1-g  environment. 
During  the  flight  mission,  three  types  of  control  strategies  will  be  used.  The 
first  will  be  the  control  derived  and  tested  on  Earth.  The  second  will  be 
derived  from  a  theoretical  0-g  model,  and  the  third  will  be  one  derived 
during  the  flight  and  uploaded  after  dynamic  tests  have  been  made  on  the 
test  article.  These  experiments  will  explore  how  accurately  closed  loop 
systems  can  be  tested  on  earth  before  being  put  into  space. 

The  control  objectives  focus  on  the  accuracy  of  the  inertial  angle  of  the 
pointing  payloads.  The  performance  objectives  are  threefold: 

*  Pointing  accuracy  of  one  payload  in  the  presence  of  disturbances. 

*  Tracking  accuracy  of  one  payload  in  the  presence  of  disturbances. 

*  Accuracy  of  both  payloads  while  both  simultaneously  point  and 
track. 
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Painting  accuracy  may  involve  a  telescope  that  needs  to  point  to  a  star 
li£^t  years  away.  One  can  easily  imagine  how  important  the  accuracy  is. 
Tracking  accuracy  may  be  necessary  on  a  structure  like  EOS,  where 
sensors  would  scan  the  earth’s  surface.  The  payload  would  then  need  to 
move  in  a  smooth  continuous  motion. 


3.3  Current  Research 

The  current  MACE  hardware  is  in  the  first  of  three  hardware  phases. 
The  Development  Model,  which  is  located  at  MPT,  is  the  first  iteration  of  the 
hardware.  The  second  iteration  will  be  called  the  Engineering  Model,  and 
the  third  iteration  is  the  flight  hardware.  There  are  currently  various 
parallel  research  projects  being  conducted  on  different  aspects  of  MACE. 
Saarmaa  [19]  has  conducted  thorough  open  loop  testing  of  the  hardware 
and  the  suspension  system.  These  tests  have  then  been  compared  to  the 
NASTRAN  finite  element  model.  Miller,  Sepe,  Rey,  Saarmaa,  and  Crawley 
[14]  have  incorporated  gravity  and  suspension  effects  into  a  finite  element 
model  of  the  development  hardware.  This  study  showed  that  low 
fundamental  frequency,  multiple  payload  devices  such  as  MACE  exhibit 
significant  suspension  and  gravity  coupling.  The  1-g  FEM  models 
displayed  much  better  correlation  to  experimental  data  than  the  0-g  models. 
Miller,  Saarmaa,  and  Jacques  [13]  have  used  measurement  based  models 
in  deriving  SISO  (single  input  single  output)  and  SITO  (single  input  two 
output)  controllers  using  Linear  Quadratic  Gaussian  design.  They  used  a 
nonlinear  curve  fitting  routine  to  derive  an  analytical  model  based  on 
experimental  frequency  response  data.  These  controllers  proved  to  be 
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effective  in  rejecting  broad  band  disturbances.  An  order  of  magnitude 
improvement  in  pointing  accuracy  was  achieved.  Other  current  research 
efforts  involve  developing  a  three  dimensional  nonlinear  model  of  MACE, 
and  the  use  of  classical  control  techniques  to  achieve  more  robustness  in 
control  design 


FIGURE  3.2:  MACE  Development  Model  Schematic. 

3.4  Development  Model  Hardware 

MACE  consists  of  a  long  flexible  bus,  a  three>axis  torque  wheel 
assembly,  and  a  two-axis  gimbal  payload.  There  are  also  various  sensors 
mounted  on  the  assembly.  The  entire  system  is  attached  to  a  sophisticated 
pneumatic  suspension  system  at  three  points  along  the  bus.  A  diagram  of 
the  hardware  and  its  components  are  shown  in  figure  3.2. 
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The  specifications  for  the  hardware  dictated  a  structure  that  is  very 
flexible  yet  durable  for  long-term  testing.  For  this  reason,  the  bus  struts  are 
made  out  of  hollow  Lexan.  There  are  four  struts  connected  by  aluminum 
collars  to  the  five  nodes.  These  collars  allow  for  easy  manual  assembling 
and  disassembling  of  the  bus.  The  nodes  are  machined  with  four  threaded 
holes  on  each  side  for  mounting  of  actuators  and  sensors.  The  bus  is 
approximately  1.5  meters  long. 

The  torque  wheel  assembly  is  used  to  control  the  attitude  of  the  bus  in 
space.  Three  inertia  wheels  are  mounted  orthogonally  on  a  triangular  base 
which  is  attached  to  node  three.  When  the  wheels  accelerate,  they  exert  a 
reaction  torque  on  the  structure.  With  precise  control  of  the  three  wheels, 
attitude  stability  in  all  three  axes  can  be  achieved.  The  unit  can  exert 
torques  of  40  oz.-in.  for  up  to  five  seconds,  and  the  Aerotech  1017  servo 
motors  can  drive  the  wheels  up  to  6000  KPM. 

A  single  gimbal  actuator  is  attached  at  node  one.  A  mass  of  steel  acts  as 
a  dummy  gimbal  on  the  other  side  of  the  bus  at  node  five  while  the 
Engineering  Model  gimbal  is  under  development.  The  two-axis  actuation  is 
achieved  through  a  yolk  design  where  one  gimbal  stage  is  contained  within 
the  other.  The  Inland  RBE  03000-B50  brushless  DC  torque  motors  directly 
drive  the  payload  in  both  axes.  These  motors  are  sized  to  exert  2000  oz.-in  of 
torque  and  to  have  a  pointing  accuracy  of  1  arcmin.  A  picture  of  the  gimbal 
without  the  payload  can  attached  is  shown  in  figure  3.3. 

There  are  also  plans  for  attaching  piezoelectric  actuators  at  strut  two  to 
provide  more  options  for  control.  This  active  strut  could  be  used  to  dampen 
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vibration  in  the  structure.  As  of  this  writing,  these  struts  were  under 
development. 


<  - 

¥ 

I 

1 


There  are  currently  28  sensors  attached  at  various  points  on  MACE. 
Sixteen  strain  gauges  (two  pairs  at  each  strut)  measure  the  bending 
deflection  in  two  directions  on  the  bus.  There  are  also  two  three-axis  rate 
gyros;  one  mounted  underneath  the  ton]ue  wheels  at  node  three,  and  the 
other  inside  the  payload  can.  These  measure  inertial  velocity  in  the  three 
directions,  and  the  signal  can  be  stably  integrated  to  acquire  the  inertial 
position.  There  are  three  single  axis  load  cells  at  the  connections  to  each 
suspension  cable,  which  measure  the  amount  of  force  being  transmitted 
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from  the  suspensiozi  system  to  the  structure.  Canon  laser  rotary  encoders 
are  attached  to  each  motor  shaft  of  the  gimbal  to  measure  the  relative  an£^e 
between  the  bus  and  the  payload.  Finally,  five  accelerometers  are  attached 
to  nodes  two  and  four  (two  at  node  two  and  three  at  node  four)  to  measure 
acceleration  in  all  three  axes.  All  of  these  sensors  are  routed  through  a 
breakout  panel  of  various  sezxsor  amplifiers  and  conditioners. 

The  structure  is  supported  by  CSA  Engineering’s  Zero-G  Suspension 
System  which  connects  to  MACE  at  nodes  1,  3,  and  5.  These  nodes  were 
chosen  to  off-load  the  three  high  mass  locations  along  the  bus.  The 
suspension  system  employs  pneumatic  pistons  and  a  feedback  loop  to 
actively  control  the  vertical  position  of  each  node.  It  allows  the  structure  to 
move  .0635m  in  the  vertical  direction  fiiction-free.  The  plunge  mode  of  the 
pneumatic  suspension  system  is  variable  and  usually  set  to  about  0.2Hz. 
The  cable  length  was  determined  such  that  the  pendulum  frequency  of  the 
suspension  system  would  not  affect  the  dynamics  of  the  plant.  Rule  of 
thumb  demands  that  this  firequency  be  less  than  an  order  of  magnitude  of 
the  first  bending  mode.  Space  limitations  constrained  the  cables  to  be  no 
longer  than  4.6m.  This  translates  to  a  pendulum  frequency  of 
approximately  0.23Hz,  reasonably  dose  to  one  tenth  of  the  predicted  bending 
mode  at  around  1.8Hz.  The  torsional  pendulum  mode  is  also  low  (*•  0.24Hz). 
There  are  also  horizontal  plane  suspension  frequences  and  violin  modes  of 
the  cables  ,  however,  which  have  higher  frequences  and  may  affect  the 
dynamics  of  the  structure.  Also,  only  the  points  where  the  suspension 
attaches  to  the  structure  are  suspended.  Other  points  in  the  structure  are 
still  in  a  gravity  field.  As  this  analysis  illustrates,  a  sophisticated 
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suspexision  system  can  be  emplojred,  but  there  are  still  numerous  problems 
in  trying  to  simulate  a  zero  gravity  environment  on  earth. 

All  of  the  control  and  time  data  acquisition  was  done  throu£^  Int^rated 
System’s  AC- 100  real  time  controller.  The  AC- 100  S3mtem  is  comprised  of  a 
DEC3100  workstation,  the  hardware  interface  and  processor  board  box,  and 
the  required  compilers  and  development  software.  The  system  has  eight 
encoder  inputs,  8  D/A  channels,  16  A/D  channels,  2  fast  analog  channels, 
and  32  parallel  digital  I/O’s.  The  main  processor  card  is  an  80386  with  a 
Weitek  3167  math  co-processor.  Integrated  Systems  has  provided  an 
interface  between  their  MatrizX  control  analysis  software  package  and  the 
processor  board  such  that  the  user  may  do  her  programming  in  the 
SystemBuild  section  of  MatrixX.  This  Amplifies  the  coding  process  in  that 
the  user  can  work  in  a  higher  level  language  (block  diagrams),  to  create  the 
download  files  to  the  processor  board.  Another  advanced  feature  of  the  AC- 
100  is  the  interactive  animation.  Users  can  directly  interface  with  the 
hardware  through  the  computer  screen.  Values  firom  the  sensors  can  be 
displayed  in  real  time  and  the  user  can  also  input  commands  finom  the 
screen  with  the  mouse.  This  removes  the  need  for  extra  monitoring  or 
input  devices.  The  user  interface  is  often  convenient  and  very  easy  to  learn 
but  not  very  powerful  in  trying  to  code  complicated  control  algorithms. 


3.5  Dynamic  Analysis 

Extensive  testing  of  the  development  model  has  been  done  to  identify  the 
dynamic  characteristics  of  the  structure.  All  of  the  data  was  taken  using  a 
Tektronix  Spectrum  Analyzer  which  fed  broad  band  white  noise  into  the 
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gimbal  motor  amplifiers.  The  sensor  readings  were  read,  amplified  and  fed 
back  to  the  spectrum  analyzer  where  a  transfer  function  was  calculated. 
This  data  was  then  analyzed  using  STAR  ID  software  to  determine  the 
resonant  frequencies  and  damping  ratios.  The  gimbal  motor  to  node  four 
vertical  accelerometer  transfer  function  is  shown  in  figures  3.4  •  3.5.  The 
fundamental  bending  frequency  is  at  1.72Hz.  Some  of  the  modes,  such  as 
the  ones  at  3.45  and  4.58Hz  were  seen  more  in  the  horizontal  plane.  This 
may  indicate  that  they  are  either  horizontal  or  torsional  modes.  A 
summary  of  the  predominant  vertical  modes  for  the  first  50Hz  is  given  in 
Table  3.1. 


TABT.E:^.!:  a  summaiy  of  the  modal  fitquencies  and  daiiq)ing  ratios 
of  the  MACS  test  article  in  the  vertical  plane. 
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Chapter  3  The  Ttstied: 


OPEN  LOOP  TRANSFER  FUNCTION 
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FIGURE  3.4;  Outer  gimbal  stage  to  node  four  accelerometer  loop 

gain  transfer  function. 


OPEN  LOOP  PHASE  TRANSFER  FUNCTION 


FIGURE  3.5:  Outer  gimbal  stage  to  node  four  accelerometer  open  lo(q> 

phase  transfer  function. 
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Chapter  3  The  ‘TtstBed: 


MACE  Development  Model 


Node  3 


Hiird  Vertical  Mode 
(8.92  Hz) 


FIGURE  3.6:  An  illustration  of  the  first  three  vertical  bending 
modes  of  the  MAdlE  develcqnnent  model 


The  mode  shapes  of  these  frequencies  can  be  predicted  by  looking  at  the 
phase  of  each  mode  and  also  by  comparing  them  to  those  derived  from  finite 
element  models.  An  illustration  of  these  mode  shapes  is  given  in  figure  3.6. 


The  frequency  identification  experiments  show  that  MACE  is  a  highly 
flexible  system  with  many  excitable  modes  with  small  damping  ratios.  The 
first  bending  mode  is  under  2Hz  and  the  three  next  bending  modes  are 
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under  lOHz.  These  dynamic  characteristics  make  MACE  an  ideal  testbed 
for  active  control  experiments,  both  closed  loop  and  open  loop  ,  i.e.  impulse 
shaping. 


3.6  Gimbal  Control 

The  gimbal  actuator  was  detached  firom  the  bus  and  mounted  on  a  test 
base  to  investigate  the  characteristics  of  the  gimbal.  The  goal  was  to 
identify  the  modal  frequencies  and  to  design  a  slewing  controller  around 
the  system.  The  specifications  for  the  original  gimbal  design  require  that 
the  gimbal  assembly  be  capable  of  achieving  a  bandwidth  of  30Hz.  There 
were  certain  challenges  in  designing  the  controller.  Some  of  the  issues  are 
dealt  with  in  this  section. 

The  outer  gimbal  stage  was  modeled  as  a  rigid  pendulum  in  a  linear 
gravity  field.  The  gimbal  is  only  capable  of  slewing  ±30®,  and  gfravity  is 
linear  to  ±4%  in  this  region.  The  gimbal  payload  arm  is  made  of  a  wide 
piece  of  aluminum  and  very  stiff  in  the  outer  stage.  Its  resonant  frequency 
is  expected  to  be  over  lOOHz  so  it  was  omitted  from  the  model.  The  dynamics 
of  the  motor  and  amplifier  were  assumed  to  be  fast  enough  to  be  modeled  as 
constant  coefRcients.  This  assumption  was  verified  by  the  vendors  specs 
and  the  successful  experiment  results.  The  parameters  for  the  pendulum 
model  of  the  gimbal  are  listed  in  table  3.2. 
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N 

^  VS:  ^ ■'  y_. 

m 

1.49  kg 

Length  from  motor 
shaft  to  of  payload 

1 

0.165  m 

Inertia  of  outer  stage 

J 

0.0441  kgm2 

Viscous  damping 

b 

0.06  N-m-s 

Motor  constant 

Kt 

0.6582  N-m/amp 

Current  amplifier 
constant 

Ki 

1  amp/volt 

TABLE  ^.2:  Parameters  used  for  gimbal  model. 


The  physical  properties  of  the  gimbal  were  measured  from  experimental 
testing.  The  motor  and  amplifier  constants  are  supplied  from  the  vendor. 
The  model  can  then  be  represented  as  a  transfer  function  in  the  form: 

flut  = - -  (3.1) 

in  Js^  +  bs+  mgl 


or  with  the  values  added: 

^  = - -M582 -  (3.2) 

V  0.0441s2  + 0.06s +  2.685 

This  simple  model  was  used  to  develop  a  PID  controller.  A  Butterworth 
pattern  of  closed  loop  poles  were  placed  and  a  bandwidth  of  3Hz  was 
indicated.  The  resulting  controller  had  gains  of:  Kp  s  22  volt/rad,  Ky  =  1.9 
volt/deg/s,  and  Ki  s  45  volt-s/deg.  The  controller  had  good  success  although 
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some  overshoot  was  seen  because  the  integrator  places  a  zero  near  the 
closed  loop  poles.  Attempts  to  increase  the  bandwidth  with  this  method 
failed  due  to  unmodeled  dynamics. 

Further  testing  revealed  that  there  were  three  different  causes  of 
instability  at  higher  gains  with  the  simple  PID  control.  The  first  was  the 
resonant  mode  of  the  structural  mount  that  holds  the  gimbal.  This  mode 
was  observed  at  around  60Hz.  It  was  identified  as  the  mount  because  the 
manual  stabilizing  of  the  mount  deleted  the  resonance  at  that  mode.  Also, 
the  addition  of  viscoelastic  material  to  parts  of  the  test  mount  resulted  in  a 
shifting  of  that  frequency.  A  second  order  notch  filter  was  designed  to 
temporarily  erase  the  problem.  A  second  unstable  mode  with  very  high 
gain  was  found  at  arotmd  460Hz.  This  mode  was  a  structural  mode  of  the 
gimbal,  most  likely  the  encoder  shaft  torsional  mode.  The  shaft  on  the  outer 
stage  is  rather  long  due  to  the  yolk  shape  of  the  design.  Since  this  is  such  a 
high  firequency,  a  low  pass  filter  was  deemed  satisfactory  for  the  problem. 
An  analog  eighth-order  Bessel  lowpass  filter  with  a  comer  frequency  of 
lOOHz  was  placed  in  the  loop  to  neutralize  the  instability.  The  Bessel  filter 
has  a  linear  phase  function,  and  therefore  causes  a  constant  time  delay  for 
all  frequencies.  This  time  delay  was  calculated  to  be  0.003  seconds. 

The  third  cause  of  instability  was  the  most  problematic.  A  dynamic 
analysis  of  the  computer  showed  that  the  processor  delay  is  on  the  order  of 
2.7  cycle  periods.  At  a  servo  rate  of  .001  seconds,  this  translates  to  a  delay  of 
.0027  seconds  for  each  feedback  signal  to  go  through  the  processor.  At  the 
desired  bandwidth  frequency  of  30Hz  this  delay  is  a  phase  loss  of  29**.  The 
phase  lag  causes  marginal  modes  to  go  unstable  very  quickly.  The  time 
delay  induces  a  linear  phase  loss  so  the  higher  firequendes  lose  more  phase 
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margin.  Along  with  the  newly  acquired  .003  seconds  of  delay  from  the 
lowpass  filters,  this  added  up  to  be  the  overwhelming  cause  of  instability  at 
high  gains.  Since  it  is  unreasonable  to  filter  out  all  the  unstable  modes,  it 
was  decided  to  include  the  computer  and  filter  dela3rs  in  the  system  model. 

A  discrete  state  space  model  was  formulated  that  indudes  the  computer 
delay  and  the  new  lowpass  filter  delay.  The  continuous  model  was  modified 
using  the  z-transform  method  to  derive  a  finite  difference  form  of  the  state 
matrix  equations.  We  start  with  the  continuous  form  of  a  system  model 
that  indudes  a  delay  in  control  action: 

x(r)  =  Fx(0  +  Gu(r  -  X)  ^2  2^ 

y  =  Hx 

where  x  is  the  state  variable,  F  is  the  dynamic  matrix,  G  is  the  input 
matrix,  H  is  the  output  matrix,  and  X  is  the  number  of  sampling  periods 
that  the  controller  delays.  As  described  in  Franklin[9],  we  can  then  derive 
the  discrete  form  of  this  model: 


x(it  +  1)  =  O  x{k)  +  Ti  u(it  -  /)  +  r2u(^  -  /  +  1) 
y{k)  =  Hx(fc) 


(3.4) 


where  it  is  an  integer  increment  in  time,  O  is  the  z-transform  of  the 
dynamic  matrix,  and  H  is  the  output  matrix.  F i  and  T2  are  wei^ted  sums 
of  the  z-transform  input  matrix  F.  The  values  of  the  matrices  are 
determined  by  the  non-integer  part  of  the  delay  X.  For  example,  if  the 
system  delay  were  a  whole  number  such  as  1  cycle  period,  Fi  would  equal  F 
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and  r2  would  equal  zero.  This  weighting  works  to  account  for  non-integer 
system  delays.  For  these  equations  to  actually  come  to  realization,  we  need 
to  introduce  I  new  variables  as  states,  one  for  each  sampling  period  of  delay. 
The  matrix  structure  of  the  equations  now  looks  like  the  following. 

■  x(k  +  1) 
x»+i(^+l) 
x^2(k+l) 

_x^i(k+l\ 


What  is  happening  is  that  each  new  state  variable  x  is  storing  a  past 
input  command  and  then  passing  that  input  command  to  the  variable  above 
it  at  every  time  step.  When  this  input  gets  to  the  “top”  of  the  matrix,  it  is 
then  multiplied  by  Fi  and  r2  to  be  used  to  evaluate  the  new  state  matrix. 
The  augmented  states  are  simply  serving  as  buffers  which  hold  past  inputs 
in  its  memory. 

The  gimbal  model  was  transformed  into  this  format  with  X=5.7  cycle 
periods  (2.7  +  3.0).  The  second  order  model  has  now  grown  to  an  eighth 
order  model.  A  full  state  controller  was  then  designed  on  this  discrete 
model.  The  eight  feedback  gains  obtained  illustrate  an  unorthodox  but 
effective  method  of  dealing  with  system  time  delays  in  control  systems.  The 
first  two  feedback  g..  .ns  are  from  the  two  original  states  of  the  system,  the 
angle  and  angular  velocity  of  the  encoder  output.  The  next  six  feedback 
gains  are  from  the  augmented  states  which  represent  the  six  past  input 
values.  The  control  law  now  takes  the  form: 


O  Fj  F2  0  **•  0 

r  xik)  1 

fol 

0  0  1  0  •••  0 

Xn+lik) 

0 

— 

0  0  0  1  ...  0 

•  *  *  •  i 

•  •  •  • 

Xn+2(k) 

+ 

0 

— 1 

0 

• 

• 

« 

0 

0 

0 

1 

LiJ 

uik) 


y(ife)  =  [H  0] 


X 

X^2 

LJr*+L 


(3.5) 


(3.6) 
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where  J  is  the  proportional  and  derivative  gains,  and  G  is  the  vector  of 
augmented  state  gains.  The  implementation  of  these  states  is  similar  to  the 
discrete  modeling  of  them.  The  calculated  inputs  are  fed  back  into  the 
control  loop  and  queued  in  another  discrete  state  space  S3rstem  for  the 
appropriate  number  of  cyde  periods  bdbre  they  are  released.  They  are 
multiplied  by  their  corresponding  gain  values  once  they  leave  the  delay 
state  space  system.  This  huffer^  ^rp®  of  control  feature  deals  directly  with 
the  phase  loss  involved  with  system  delays.  A  blodc  diagram  of  the  control 
algorithm  is  shown  in  figure  3.7.  The  inputs  are  the  encoder  counter  and 
three  command  inputs  fi:om  the  user  which  indicate  the  two  angles  that 
she  wants  to  dew  between.  The  output  from  the  controller  is  the  command 
voltage  into  the  gimbal  motor  amplifier.  This  voltage  is  cut  off  at  ±10  volts  to 
avoid  saturation. 

A  bandwidth  of  15Hz  was  obtained  using  this  strategy.  Subsequent 
experimentation  with  notch  filters  obtained  the  30Hz  specification 
bandwidth.  The  transfer  functions  for  the  15Hz  bandwidth  controller  with 
the  buffer  delay  function  are  shown  in  figures  3.8  -  3.9.  These  are  compared 
to  the  theoretical  results  predicted  firom  the  discrete  augmented  state 
model.  Notice  that  the  unmodeled  modes  have  altered  the  bandwidth 
somewhat.  There  is  still  noticeable  phase  loss  due  to  delays,  but  the  delay 
states  have  yielded  an  improvement  over  the  previous  phase  plot.  The 
implementation  of  the  buffers  had  a  realizable  stabilizing  effect  on  the 
controller.  A  summary  of  the  gains  for  the  15Hz  bandwidth  controller  is 
shown  in  table  3.2. 


MatrixX  SystemBuild  Diagram 
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FIGURES.?:  A  block  diagram  of  control  algorithm  for  a  PD 
contioller  with  augmented  delay  buffer  states. 
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CLOSED  LOOP  GAIN  UtANSFER  FUNCnON 


FIGURE  3.8:  Qosedloq)  gain  txansferfiinctira  of  command  to 

encoder  angle. 


CLOSED  LOOP  PHASE  IRANSFER  FUNCTION 
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552  volts/rad 

Kv 

11.6  volts/rad/sec 

gl 

0.180  volts/volt 

g2 

0.249  volts/volt 

g3 

0.237  voltsA^olt 

g4 

0.224  voltsArolt 

gs 

0.212  voltsArolt 

g6 

0.200  volts/volt 

TART.F.  A  summary  of  the  parameters  for  ISHz  bandwidth  COTtroller 

with  delay  buffer  states. 


3.7  Difficulties  Encountered  and  other 
Considerations 

Dealing  with  computer  and  filter  delays  was  a  good  exercise  in  control 
theory  but  an  uimecessary  one.  Further  research  on  the  MACE  test  article 
will  demand  much  more  computing  power  and  will  tolerate  less  delay  than 
these  initial  experiments.  Higher  order  controllers  will  require  more 
computation  and  may  decrease  the  servo  rate.  Since  the  processor  delay  is 
a  function  of  the  servo  rate,  it  will  become  increasingly  more  troublesome. 
Plans  for  a  new  processor  board  may  help  the  situation.  Also,  the  instabiP 
at  460Hz  is  a  design  oversight  that  needs  to  be  resolved.  A  solution  may 
erase  the  need  for  a  filter  in  the  loop  and  cause  less  phase  loss.  This  may 
come  about  in  the  second  generation  Engineering  Model  gimbal  which  has 
taken  this  into  consideration. 
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There  are  some  other  control  issues  for  the  gimbal  that  were  not  covered. 
An  integrator  feedback  was  not  used  because  a  reasonable  method  of 
modeling  the  integrator  in  the  augmented  state  model  was  not  found. 
However,  ad  hoc  methods  of  prediction  and  trial  and  error  may  prove  to  be 
an  effective  way  to  add  an  integrator  term  to  the  control  law.  Another 
method  of  reducing  the  steady  state  error  is  to  use  a  gravity  compensator. 
This  was  also  examined  but  neglected  due  to  the  problems  with  system 
delays.  Two  methods  of  implementation  are  possible.  A  feedback  method 
would  read  the  encoder  angle  and  transmit  the  right  amount  of  current 
into  the  motor  to  hold  the  gimbal  steady  at  that  angle.  A  feedforward 
method  would  translate  the  input  command  into  a  voltage  corresponding  to 
that  angle  and  feed  it  directly  into  the  gimbal.  Both  of  these  methods  would 
increase  pointing  accuracy  and  would  work  in  addition  to  the  controller  in 
place.  A  more  legitimate  method  of  control  derivation  may  be  to  model  the 
unstable  structural  modes  of  the  gimbal  and  the  testing  mount,  rather  than 
filtering  them.  This  would  result  in  a  more  complicated  model  and  a 
higher  order  controller  that  may  further  tax  the  computer.  The  controllers 
implemented  in  this  section  were  found  to  be  sufficient  for  the  required 
slewing  and  pointing  applications. 


3.8  Closure 

In  this  chapter  we  have  covered  the  Middeck  Active  Control  Experiment 
as  the  testbed  for  the  experiments  presented  in  this  paper.  MACE  is  a 
NASA  funded  project  whose  purpose  is  to  develop  the  hardware  and  science 
specifications  for  experiments  on  board  the  Space  Shuttle  during  the 
summer  of  1994.  The  goal  of  the  project  is  to  investigate  issues  concerning 
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active  control  on  a  flexible  structure  in  a  zero  gravity  environment.  The 
Development  Model  (DM)  hardware  located  at  MTT  has  been  discussed  in 
detail.  The  suspension  system  was  seen  as  an  advanced  tool  for  simulating 
a  zero  gravity  environment,  but  not  sufficient  as  to  be  able  to  ignore  gravity 
effects.  Some  results  of  open  loop  testing  were  given  to  describe  the  dynamic 
characteristics  of  the  DM.  Finally,  a  detailed  description  of  the  procedure 
for  deriving  the  controller  for  the  gimbal  actuator  was  presented. 

Impulse  shaping  is  one  of  many  different  parallel  research  efforts  going 
on  at  the  Space  Engineering  Research  Center.  MACE  will  provide  a  vehicle 
for  the  impulse  shaping  technology  to  gain  credibility  by  securing  real  on- 
orbit  test  results.  The  goals  and  specifications  of  the  MACE  project  are 
congruent  with  those  of  my  research.  It  is  a  most  suitable  platform  for 
advanced  experiments  on  impulse  shaping. 


4  l^Smutatufn:  DISCOS 


THE  SIMULATION:  DISCOS 

_ Chapter  4 


4.1  Introduction 

There  were  two  reasons  we  wanted  to  utilize  simulation  as  a  tool  in  this 
research.  The  first  is  that  impulse  shaping  had  been  verified  in  simple 
S3rstems  and  linear  models,  but  nonlinearities  had  not  been  addressed.  We 
wanted  to  present  a  comparison  between  the  linear  and  nonlinear 
simulations.  The  hope  was  to  gain  some  insight  into  how  impulse  shaping 
afiects  nonlinear  systems.  Second,  simulation  was  a  method  of  verifying  our 
algorithms  before  applying  them  to  the  hardware.  Nonlinear  simulation  was 
a  tool  that  was  available  before  the  hardware.  This  chapter  will  briefly 
describe  the  DISCOS  simulation  code  and  the  corresponding  model  of  MACE. 
We  will  end  the  chapter  with  some  results  of  dynamic  analysis  and  impulse 
shaping  results.  Bodley  [3]  should  be  consulted  for  a  more  in  depth  review  of 
DISCOS. 
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DISCOS  stands  for  Dynamic  Interaction  Simulation  of  Controls  and 
Structure.  It  is  a  fairly  large  program  written  in  FORTRAN  developed  by 
Martin  Marietta  and  Goddard  Flight  Center  in  the  1970’s.  Its  purpose  is  to 
provide  a  tool  for  studying  the  nonlinear  behavior  of  large,  miiltiple-body, 
flexible  space  structures.  It  formulates  and  calculates  time  domain  solutions 
to  the  nonlinear  differential  equations  of  motion.  It  also  has  a  linearization 
capability  for  frequency  domain  and  control  stability  analysis.  Some  features 
which  allow  it  specifically  to  be  used  in  space  systems  are  the  inclusion  of 
gravity  gradients,  thermal  gradients,  solar  pressure,  and  drag.  It  is  an  ideal 
program  for  analyzing  large  satellites  or  space  vehicles.  We  decided  to  use  it 
to  simulate  MACE. 

The  MACE  model  input  into  DISCOS  is  derived  from  an  early  version  of 
an  ADINA  finite  element  model  developed  at  SERC.  The  details  of  the 
derivation  of  this  model  are  described  in  Padilla  [17].  DISCOS  sees  a  system 
as  a  collection  of  bodies  connected  by  hinges.  The  bodies  may  be  flexible  or 
rigid  and  the  hinges  can  be  defined  as  clamps,  pins,  rollers,  or  spring  damper 
systems.  MACE  was  modeled  as  three  bodies,  two  rigid  payloads  and  a 
flexible  bus  coimected  by  pins.  A  schematic  of  the  MACE  model  is  shown  in 
figure  4.1. 

The  bus,  two  gimbals,  and  torque  wheels  are  modeled  as  one  long  beam 
with  concentrated  masses  at  the  actuator  locations.  At  the  center  of  the  bus 
is  an  actuator  point  where  outside  torques  and  forces  can  be  applied.  This  is 
how  we’ll  simulate  the  action  of  the  torque  wheels.  The  gimbal  actions  will  be 
simulated  by  a  hinge  torque  between  the  bus  and  the  payloads.  A  sensor  is 
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placed  at  the  hinge  between  payload  two  and  the  bus  to  provide  a 
measurement  of  vibration. 


MACE  DISCOS  model  schematic 

4.55  kg 


Sensor 


Body  1  =  Bus  -f  Gimbals  +  Torque  Wheels 
Body  2  =  Payload  1  (Rigid) 

Body  3  =  Payload  2  (Rigid) 


FIGURE  4.1:  A  schematic  of  what  MACE  looks  like  in  DISCOS. 


The  flexibility  in  the  bus  is  modeled  as  a  set  of  elastic  displacement 
coordinates  projected  onto  the  inertial  axis  which  are  derived  from  the 
vibration  modal  information  from  the  FEM  model.  These  modes,  along  with 
the  mass  and  inertia  information  for  each  point  on  the  bus,  are  defined 
kinematically  in  DISCOS  with  respect  to  the  body-fixed  reference  frame  on 
the  bus.  The  first  nine  modes  of  the  system  were  modeled. 
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4.3  Dynamic  Analysis 

To  find  the  modal  firequendes  of  the  full  model  in  different  configurations, 
frequency  identification  tests  were  performed.  Wide  band  white  noise  was 
generated  fiinm  the  computer  and  used  as  input  into  the  payload  one  gimbal. 
In  DISCOS,  this  was  translated  as  a  hinge  torque  between  bodies  one  and 
two.  The  sensor  output  from  the  sensor  between  the  bus  and  payload  two 
was  read.  Ten  simulations  were  completed  and  an  average  of  the  Fast 
Fourier  Transform  of  the  runs  was  analyzed  to  identify  the  resonant  modes. 
Several  configurations  were  tested  to  distinguish  the  differences  in  modal 
firequendes.  PD  controllers  were  added  both  payloads  and  at  the  attitude 
controller  to  provide  a  dosed  loop  analysis  as  well.  A  summary  of  the  modal 
information  for  different  configurations  is  given  in  figure  4.2. 


KEY:  O  CONTROLLER  A  PAYLOAD  □  LOCKED 

FIGURE  4.2:  DISCOS  model  modal  firequendes  in  four 
different  configuraticnis. 
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The  parameters  used  in  this  model  come  from  an  earlier  design  of  MACE. 
Therefore,  the  dynamic  characteristics  may  be  slightly  different  from  the 
actual  hardware.  The  closed  loop  frequencies  are  also  different  because  the 
controllers  are  altered  along  with  the  system.  The  first  three  open  loop 
bending  frequencies  of  the  DISCOS  model  were  1.65Hz,  5.4Hz,  and  6.9Hz. 
The  real  modal  frequencies  were  1.72Hz,  6.70Hz,  and  8.92Hz. 

The  kinematic  nonlinearities  in  the  DISCOS  equations  allow  the 
identification  of  the  shifts  in  the  modal  fi^quendes.  These  shifts  in  firequency 
were  a  focus  in  the  analysis  of  the  simulations.  We  can  see  how  the  first  four 
frequendes  change  with  the  different  configurations.  It  was  not  clear  how 
effective  impulse  shaping  woiild  be  as  the  firequendes  change  during  a  large 
angle  slew.  Also,  verification  of  the  impulse  shaping  algorithm  was  sought  in 
a  nonlinear  formulation  of  the  equations  of  motion.  An  impulse  shaper  only 
seeks  the  modal  information  from  a  system  (cd  and  Q.  Since  this  assumption 
comes  from  fundamental  linear  theory,  it  was  desired  to  see  its  effectiveness 
on  a  nonlinear  system. 


4.4  Simulation  Results 

The  first  set  of  tests  completed  were  open  loop.  The  payloads  were 
damped  in  the  vertical  position  to  the  bus.  The  attitude  controller  was  used 
as  the  actuator  of  the  torques  on  the  bus.  A  bang-bang  tyi>e  of  input  was  used 
to  torque  the  system  to  an  angle  and  stop  it.  The  inertial  angle  of  the  bus 
was  used  as  an  indicator  of  the  system  vibration.  The  input  was  then  shaped 
tising  a  three  impulse  sequence  for  the  first  mode  at  1.65Hz.  The  first  mode 
dominated  the  response  of  the  system  and  the  suppression  of  that  frequency 


ANGLE  (deg)  TOR<^{Nm) 


TIME(mc) 

FIGURE  4.3:  Shaped  and  unshipped  tmque  inputs  into  atdn^ 
c^mtroller.  1  mode  (1.65  Hz) 


BUS  ANGLE  RESPC»<SE  TO  UNSHAFED  TORQUE 


TIME(aec) 


FIGURE  4.4:  Bus  inertial  angle  response  to  unshaped  torque 
input  into  attitude  controller. 
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•1  0  1  2  3  4  5 

TIME(mc) 

FIGURE  4.5;  Bus  inerdal  angle  response  to  shaped  torque 
input  into  attitude  controUer  (l.tiSHz) 

would  cancel  most  of  the  vibration.  The  same  output  was  recorded  and 
compared  to  the  unshaped  case.  A  sunmiary  of  the  results  is  given  in  figure 
4.3  -  4.6. 

The  unshai)ed  response  shows  the  multiple  mode  behavior  of  the  system. 
It  is  clear  that  the  bang-bang  input  has  excited  much  vibration  in  the  bus. 
The  shaped  response  shows  a  vast  improvement  over  the  unshaped  response. 
There  is  almost  no  visible  residual  vibration  in  the  first  mode.  Notice  how 
the  input  is  now  staircased  for  the  three  impulses  in  the  sequence.  The  time 
delay  incurred  firom  the  shaping  is  about  0.606  seconds.  The  settling  time  of 
the  vibration  however  is  almost  zero. 
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A  dosed  loop  simulatioii  was  the  goal  of  the  next  test.  PD  controllers  at 
the  three  actuators  provided  a  configuration  that  is  more  like  what  will 
happen  on  the  test  flight.  Payload  one  was  given  a  position  command  that 
moved  fix)m  -20**  to  +20**.  The  deflection  in  the  Y  direction  was  read  from  the 
sensor  mounted  at  payload  two.  The  command  was  then  shaped  for  a  direct 
solution  sequence  for  the  first  three  modes.  A  summary  of  the  results  is 
given  in  figures  4.7-4.8.  Notice  that  the  sensor  reads  motion  relative  to  the 
body-fixed  reference  fhune  of  the  bus  and  hence  the  steady  state  deflection  is 
zero.  The  deflections  are  small  due  to  the  actions  of  the  control  system. 


PAYLOAD  ANGLE  RESPONSE  TO  SLEW 


FIGURE  4.6:  Payload  angle  response  to  40**  shaped  and  unshaped 
slews.  3  modes  (1.85, 7.5, 15.5  Fb). 
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ENDPOINT  DEFLECmON  TO  UNSHAPED  SLEW 


TIME  (sec) 


FIGURE  4.8:  Sensor  deflection  in  y  direction  for  shaped 
40®  slew.  3  modes  (1.85, 7.5, 15.5  Hz) 
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The  results  are  similar  to  that  of  the  open  loop  tests.  The  shaped  response 
shows  a  dear  improvement  in  vibration  response.  However,  a  doser  look  at 
the  responses  will  reveal  that  the  vibration  response  is  not  absolutely  zero. 
We  notice  that  in  the  slew  from  -20  to  -f  20  degrees  the  first  modal  firequency 
has  shifted  from  1.80Hz  to  2.01Hz.  The  shaper  has  been  designed  for  the 
vertical  configuration  with  the  first  mode  at  1.85Hz.  The  changes  in 
frequency  due  to  the  kinematic  nonlinearities  are  causing  a  slight  decrease  in 
the  effectiveness  of  the  impulse  shaper. 

Although  the  shaper  did  not  meet  the  absolute  zero  vibration  constraints 
set  by  the  original  equations,  it  did  show  a  significant  improvement  in  the 
vibration  response  over  unshaped  moves.  It  wasn’t  expected  that  the 
nonlinear  model  would  conform  completely  to  impulse  shaping.  More 
analysis  of  the  nonlinearities  will  be  discussed  in  the  hardware  experiments. 
These  simulation  experiments  have  proven  that  impulse  shaping  can  be 
effective  on  nonlinear  systems.  This  will  serve  as  a  basis  for  the  results  from 
hardware  tests  presented  in  the  next  chapter. 


EXPEREVIENTAL  RESULTS 

Chapter  5 


5.1  Introduction 

Data  from  the  MACE  hardware  was  sought  in  order  to  verify  the 
impulse  shaping  effects  on  a  real  system.  We  also  wanted  to  address  the 
issues  that  would  arise  from  implementation.  The  experiments  presented 
in  this  chapter  are  from  the  MACE  Development  Model  and  its  supporting 
electronics  and  suspension  system.  This  chapter  will  also  provide  an 
analysis  of  the  data,  its  implications,  and  how  to  deal  with  some  of  the 
problems. 

5.2  Defining  Input  Shaping  Effectiveness 

Before  we  can  start  evaluating  experimental  results,  we  need  to  define 
what  we’re  looking  for.  There  needs  to  be  a  performance  metric  to  gauge 
the  effectiveness  of  a  shaper  in  reducing  vibration.  This  performance 
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metric,  of  cotirse,  depends  on  the  particular  application  and  is  a  function  of, 
among  other  variables,  the  percentage  of  vibration  reduction. 

The  MACE  flight  plan  calls  for  the  pointing  of  one  payload  while  the 
other  slews  to  a  different  position.  The  stability  of  the  first  payload  is  crucial 
since  the  payload  is  likely  to  be  a  precise  telescope  or  laser  pointing  at  a  far 
away  target.  The  primary  goal  then  is  to  design  a  regulator  that  rejects  all 
the  disturbances  caused  by  the  slew.  We  are  not  concerned  with  that  here. 
The  secondary  goal  is  to  be  able  to  slew  the  second  payload  without  causing 
much  disturbance  to  the  other  end  of  the  bus.  This  is  where  impulse 
shaping  is  employed. 

The  appropriate  sensor  needs  to  be  chosen  to  evaluate  the  relevant 
response  of  the  system.  The  vertical  accelerometer  at  node  four  on  the  other 
end  of  the  bus  (figure  3.2)  was  chosen  as  it  relates  closely  to  the  amount  of 
force  that  is  exerted  on  the  other  payload.  It  also  shows  good  amplitude  for 
the  mode  shapes  of  those  modes  under  lOHz.  Notice  that  unless  a  sensor  is 
mounted  on  a  node  point  of  a  particular  mode  shape,  the  same  resonances 
should  be  seen  at  any  sensor  on  the  structure.  There  may  be  additional 
zeroes  in  the  transfer  function  to  some  sensors,  but  the  poles  do  not  change. 
Singer  has  shown  that  numerator  dynamics  do  not  matter  in  impulse 
shaping.  So,  the  choice  of  sensors  just  allows  us  to  more  dosely  examine  a 
certain  point  on  the  structure;  it  does  not  determine  the  validity  of  the 
results.  Other  sensors  were  observed  to  verify  the  accuracy  of  the  first. 

MACE  was  designed  to  have  many  lightly  damped  modes  imder  20Hz  as 
seen  in  section  3.5.  In  vibration  control,  we  are  mostly  concerned  with  low 
frequency  modes.  Generally  speaking,  the  slow  modes  tend  to  last  longer 
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and  have  higher  amplitudes,  and  hence,  cause  more  of  a  delay  by 
increasing  the  settling  time.  We  arbitrarily  decided  to  concentrate  on  those 
modes  under  lOHz,  ignoring  the  higher  modes  as  less  significant.  In  a 
typical  slew,  these  modes  were  the  most  visible  and  the  cause  of  most  of  the 
vibration  problems.  In  the  vertical  plane,  6  modes  were  identified  under 
lOHz,  and  it  was  felt  that  was  a  suitably  complex  problem  for  the  impulse 
shaping  calculation  methods  available. 

The  performance  metric  becomes,  then,  the  amount  of  vibration  a 
certain  shaper  can  suppress  in  the  vertical  acceleration  at  the  other  end  of 
the  bus  for  the  modes  within  the  first  lOHz  of  the  MACE  test  article.  As  a 
general  rule,  we  will  describe  "good”  performance  as  5%  residual  vibration. 
We  then  weight  the  amount  of  vibration  suppression  against  the  complexity 
of  the  impulse  sequence  and  the  length  of  time  delay  caused  by  the 
sequence.  The  time  delay  is  a  cost  since  a  very  slow  move  regardless  of 
shape  will  always  cause  little  vibration.  The  solution  complexity  is  still  an 
issue  because  the  formulation  of  the  sequences  has  not  been  fiilly  resolved. 


5.3  Open  Loop  Results  and  Discussion 

The  first  set  of  tests  done  on  MACE  were  open  loop.  An  entered  user 
command  was  converted  to  a  current  command  that  was  directly  input  to 
the  motor.  Theoretically,  the  controller  state  of  the  plant  does  not  matter  in 
impulse  shaping.  As  seen  in  the  DISCOS  simulations  in  Chapter  4, 
shapers  can  be  just  as  effective  for  open  loop  slews  as  for  closed  loop  slews. 
The  open  loop  shaping  results  in  actual  hardware,  however,  were  very  poor 
and  showed  little  improvement  in  residual  vibration.  This  was  due  to  the 
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hardware  inconsistencies  which  were  too  severe  to  overcome  with  an  open 
loop  control  scheme.  The  slew  profiles  and  the  plant  responses  were  both 
very  inconsistent  firom  run  to  run.  An  impulse  shaper  requires  very 
precise  movement  of  the  actuator  to  be  effective.  A  shaped  command  tends 
to  be  more  irregular  than  the  original  commanded  move  and  thus  requires 
more  precision  firom  the  actuator.  There  are  many  factors  in  the  setup  that 
hinder  this  accurate  movement.  There  is  a  lot  of  steady  state  error  in  the 
gimbal  movement.  The  brushless  motor  also  has  problems  doing  consistent 
moves  for  a  given  input  current.  There  are  also  other  factors  such  as 
friction,  cogging  torque,  amplifier  dynamics,  etc.  A  more  precise  actuator 
should  produce  just  as  consistent  results  as  the  closed  loop  experiments 
which  are  covered  next. 


5.4  Closed  Loop  Frequency  Response 

5.4.1  Experimental  setup 

A  lOHz  bandwidth  PD  loop  similar  to  the  one  that  was  designed  in 
Chapter  3  was  implemented  on  the  gimbal  as  mounted  on  the  bus.  It 
should  be  noted  that  this  is  no  longer  a  full  state  feedback  compensator  and 
there  are  many  unmodeled  modes  of  the  bus  within  the  bandwidth  of  the 
controller  that  are  not  being  controlled.  The  compensator  was  found  to  be 
stable  and  robust,  however,  and  adequate  for  the  slewing  application. 

This  loop  was  closed  and  a  transfer  function  was  then  taken  from  the 
command  gimbal  input  to  the  vertical  acceleration  at  node  four.  The  setup 
for  this  experiment  is  shown  in  figure  5.1.  White  noise  is  generated  from  a 
Fourier  analyzer  and  input  into  the  real-time  controller.  This  voltage  is 


AC-100  real-time 
controller. 


FIGURE  5.1:  Diagram  of  Experimental  Setup 
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translated  to  angle  conunands  into  the  control  loop.  This  then  drives  the 
gimbal,  causing  a  response  in  the  bus  and  the  node  four  vertical 
accelerometer.  The  accelerometer  signal  is  amplified  and  input  into  the 
Fourier  analyzer.  The  anal3rzer  then  calculates  ten  averages  of  the  Fast 
Fourier  Transform  of  the  accelerometer  output  divided  by  the  transform  of 
the  white  noise,  resulting  in  a  closed  loop  frequency  response  of  the  plant. 

This  loop  was  closed  and  a  transfer  function  was  then  taken  from  the 
command  gimbal  input  to  the  vertical  acceleration  at  node  four.  The  setup 
for  this  experiment  is  shown  in  figure  5.1.  White  noise  is  generated  from  a 
Fourier  analyzer  and  input  into  the  real-time  controller.  This  voltage  is 
translated  to  angle  commands  into  the  control  loop.  This  then  drives  the 
gimbal,  causing  a  response  in  the  bus  and  the  node  four  vertical 
accelerometer.  The  accelerometer  signal  is  amplified  and  input  into  the 
Fourier  analyzer.  The  analyzer  then  calculates  ten  averages  of  the  Fast 
Fourier  Transform  of  the  accelerometer  output  divided  by  the  transform  of 
the  white  noise,  resulting  in  a  closed  loop  firequency  response  of  the  plant. 


The  results  are  shown  in  figures  5.2  -  5.3.  We  note  the  similarity  to  the 
open  loop  transfer  functions  shown  in  Chapter  3.  The  amplitudes  are 
different  because  the  units  have  changed  from  m/s^/volt  to  m/s^/deg.  The 
peaks,  however,  retain  the  same  sliapes.  This  is  due  to  the  large  relative 
size  of  the  bus  compared  to  the  gimbal  actuator.  We  can  imagine  that 
adding  a  controller  is  similar  to  adding  a  small  spring  and  damper  at  the 
first  node.  This  new  stiffiiess  has  caused  certain  modes  to  change.  Two  of 
the  visible  alterations  are  the  smaller  amplitude  and  sharper  peak  of  the 
first  mode  and  the  larger  amplitude  of  the  fifth.  The  second  vertical  mode 
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E5.2:  Qosed  loop  gain  transfer  functitm  fiiom  commanded 
gimbal  angle  to  vertical  acceleration  at  node  4. 
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at  5.23Hz  has  been  most  affected  by  the  controller  as  it  has  shifted  from  its 
prior  6.82Hz. 

5.4JS  Shaped  response 

The  experiment  was  then  repeated  with  an  impulse  shaper  placed  in 
between  the  Fourier  analyzer  and  the  controller  to  examine  the  shaping 
effects  on  the  system  in  the  frequency  domain.  The  shaper  is  programmed 
into  the  microprocessor  and  designed  to  convolve  any  input  with  the 
impulse  sequence  in  real  time.  The  white  noise  is  now  shaped  before  it 
enters  the  control  loop  and  the  shaper  becomes  a  part  of  the  plant  whose 
transfer  function  is  being  taken.  The  shaped  closed  loop  transfer  function 
is  then  acqtiired  as  shown  in  figures  5.4  -  5.5.  Results  for  two  different 
shapers  are  shown:  a  three  impulse  shaper  for  the  first  mode  at  1.75Hz  and 
a  nine  impulse  direct  solution  sequence  for  the  first,  fourth,  fifth,  and  sixth 
modes  (1.75, 5.23, 8.68,  and  9.13Hz).  These  shapers  were  chosen  to  suppress 
vibration  at  the  slowest  mode,  and  then  the  next  three  with  the  highest 
amplitudes.  The  second  and  third  modes  were  omitted  because  they  are 
primarily  horizontal  modes  as  shown  in  section  3.5. 

We  can  immediately  see  what  the  shaper  is  doing.  The  three  impulse 
sequence  for  the  first  mode  has  cut  the  amplitude  of  vibration  in  uiat  mode 
by  almost  two  orders  of  magnitude.  Moreover,  we  can  see  that  the  areas 
around  that  mode  are  also  suppressed  as  well.  This  robustness  is  even 
more  evident  in  the  direct  solution  sequence.  We  see  that  most  of  the  area 
under  the  curve  between  around  5Hz  and  13Hz  has  been  destroyed.  We  also 
notice  from  these  plots  the  repeating  nature  of  these  shapers.  They  happen 
to  cancel  modes  which  are  multiples  of  the  ones  being  shaped.  For 
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SHAPED  CLOSED  LOOP  TRANSFER  FUNCTION  (1  MODE) 


- SHAPED 

- UNSHAPED 


FREQUENCY  (HZ) 

Closed  loop  transfer  function  with  1  mode,  3  inq>ulse 
shaper  (USHz) 


SHAPED  CLOSED  LOOP  TTIANSFER  FUNCTION  (4  MODES) 
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Qosed  loop  transfer  function  with  4  mode  direct  solution 
sequence  (1.75,  5,23,  8.68,  913Hz). 
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example,  in  the  three  impulse  one  mode  shaper,  by  coincidence  we  have 
also  canceled  most  of  the  fourth,  fifth,  and  sixth  modes  as  well.  It  is  not  yet 
quite  understood  where  the  direct  solution  sequences  will  repeat,  but  it  is 
clear  that  they  have  a  unique  pattern  of  repetition.  Although  we  don’t 
cancel  as  much  amplitude  or  achieve  as  much  robustness  in  higher  modes 
by  using  technique,  we  will  see  later  that  it  can  be  just  as  effective  as 
independently  shaping  for  those  modes. 


5.4.3  Experimental  insensitivity  curves 

Recall  how  we  calculate  an  insensitivity  curve.  It  is  a  measure  of  the 
amplitude  of  vibration  response,  at  each  frequency  on  the  x  axis,  to  a  given 
shaped  move  as  a  percentage  of  the  response  to  an  unshaped  move.  An 
amplitude  of  one  would  mean  that  the  shaped  move  causes  the  same 
amount  of  residual  vibration  at  that  fmquen(7  as  the  unshaped  move.  The 
transfer  functions  we  have  taken  have  enough  data  on  them  to  produce 
experimental  versions  of  these  insensitivity  curves.  By  dividing  at  each 
point  in  frequency  the  amplitude  of  the  shaped  response  by  the  amplitude  of 
the  unshaped  response,  we  get  a  measure  of  the  percentage  of  vibration  that 
the  shaped  excitation  has  left  relative  to  the  imshaped  excitation.  The 
resulting  curve  is  in  effect  a  plot  of  the  amount  of  vibration  a  shaper  can 
cancel  for  a  range  of  frequencies.  This  can  be  directly  correlated  to  the 
theoretical  insensitivity  curves  that  have  been  produced  using  the  impulse 
response  equations.  The  experimental  insensitivity  curves  for  both  the  one 
mode  three  impulse  sequence  and  the  four  mode  direct  solution  sequences 
are  shown  in  figures  5.6  -  5.7. 


ferimntaC  ^^uCts 


EXPERIMENTAL  INSENSmYITY  CURVE  (1  MODE) 


1  mode  (1.75Hz). 


EXPERIMENTAL  INSENSITIVITY  CURVE  (4  MODES) 


LZ:  Insensitivity  curve  for  4  mode  direct  solutitm 
sequence  (1.75,  5.23,  8.68,  9.13Hz). 
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These  curves  provide  a  more  visual  measure  of  the  robustness  of  these 
shapers  to  errors  in  plant  parameter  estimation.  For  the  first  mode,  we  can 
see  we’re  within  5%  vibration  error  in  the  range  of  1.5Hz  to  2.0Hz.  We  also 
see  the  repetition  of  the  zero  vibration  constraint  at  odd  multiples  of  l.TSHz. 
The  loss  of  effectiveness  at  higher  modes  is  a  result  of  the  digitization  of  the 
impulse  sequence  (we  can  only  approximate  the  times  of  the  impulses  as 
acctirately  the  servo  rate  of  the  controller  allows;  in  this  case  .001  sec).  The 
effect  of  digitization  is  seen  more  in  the  higher  repeated  modes  as  can  be 
seen  in  the  figures.  In  the  four  mode  direct  solution  plot,  we  see  the  same 
robustness  in  the  1.75Hz  range,  and  we  see  even  more  robustness  in  the  8  to 
lOHz  range  where  we  have  shaped  for  two  modes.  This  is  a  case  where  the 
robustness  of  the  shaper  for  two  modes  which  are  dose  together  “overlap” 
and  cause  a  large  region  to  be  very  robust.  That  whole  region  appears  to 
have  close  to  zero  percent  vibration.  We  also  note  that  the  robustness  of  the 
5.23Hz  region  is  about  the  same  as  that  of  the  1.75Hz  region.  As  discussed 
in  Chapter  2,  this  is  due  to  the  direct  solution  formulation.  A  convolved 
sequence  would  show  approximately  three  times  more  robustness  in  the 
5.23Hz  range.  For  the  time  savings  gained  fi:t)m  the  direct  solution  method, 
however,  this  sequence  seems  to  be  a  more  efficient  one. 


When  the  experimental  results  are  compared  to  the  theoretical  data,  we 
see  a  very  close  correlation.  Equation  2.3  is  used  to  calculate  at  each  point  in 
frequency  the  amplitude  of  vibration  to  obtain  the  theoretical  insensitivity 
curves  for  the  particular  sequences.  Both  sets  of  data  are  shown  in  fig^ures 
5.8-5.9.  The  theoretical  data  has  also  been  calculated  using  the  digitized 
sequences.  We  see  that  the  experimental  data  not  only  follows  the  trends  of 
the  theoretical  data,  but  even  mimics  the  small  details,  such  as  the  tiny 
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EXPERIMENTAL  VS.  THEORETICAL  INSENSITIVITY 


0  2  4  6  8  10  12  14  16  18  20 


FREQUENCY  (HZ) 

FIGUi-  “  5.8:  A  comparison  of  experimental  versus  theoretical 
insensitivity  curves  for  1  nxxie  3  impulse  shaper  (1.7SHz) 


EXPERIMENTAL  VS.  THEORETICAL  INSENSITIVITY 
CURVES  (4  MODES) 


FIGURE  5.9:  A  comparison  of  experimental  versus  theoretical 
insensitivity  curves  for  4  mode  direct  solutitxi  sequence 
(1.75,  5.23  8.68,  9.13H2). 


mtaC^gsults 


88 


hump  at  12.5Hz  for  the  four  mode  direct  solution  sequence.  The 
experimental  data  obtained  through  this  procedure  indicates  a  good 
correlation  between  the  theoretical  model  and  the  actual  hardware. 


5.4.4  Why  does  it  work  so  well? 

Any  time  experimental  data  looks  this  good,  one  must  question  why  it 
does  so.  The  theoretical  data  is  based  on  a  linear  representation  of  the 
plant.  It  is  derived  from  the  second-order  linear  response  of  a  system  to  a 
series  of  impulses.  We  can  deduce  from  the  experimental  data  that  the 
system  behaves  like  the  linear  system  described,  or  a  dose  approximation  of 
one.  There  are  a  couple  of  reasons  that  MACE  may  look  linear  in  this 
configuration. 

When  the  frequency  data  is  being  taken,  the  plant  does  not  move  much. 
The  white  noise  amplitude  is  kept  small  due  to  the  amplifier  constraints 
and  there  is  no  DC  offset,  so  there  is  no  rigid  body  movement.  There  is 
enough  movement  to  induce  vibration  in  the  bus,  and  enough  amplitude  to 
see  high  coherence  in  the  frequency  response  data.  However,  the  envelope 
of  vibration  in  the  bus  is  small  and  the  gimbal  angle  displacement  is 
minimal  (<  Typically,  Dynamic  systems  that  have  small  angles  or 

displacements  can  be  modeled  as  linear  systems.  So  the  small  amplitudes 
of  the  test  lend  itself  to  exhibiting  very  linear  behavior. 


Another  reason  for  the  good  correlation  is  the  lightly  damped  modes  of 
MACE.  When  the  impulse  sequences  and  the  insensitivity  curves  were 
calculated,  we  assumed  that  every  mode  had  zero  percent  damping.  This  is 
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because,  in  order  to  calculate  the  response  of  a  lErequency  to  the  impulse 
sequence,  we  need  to  know  the  damping  ratios  at  each  frequency.  This  is 
difficult  to  know,  so  we  generally  assume  that  th>  imping  ratios  at  each 
frequency  are  the  same  as  that  of  the  most  lightly  damped  mode  being 
shaped.  The  modes  within  the  first  20Hz  of  MACE  are  veiy  lightly  damped; 
all  under  2%  damping.  So  the  zero  percent  damping  was  a  good 
assumption  in  this  case.  Each  firequency  imder  20Hz  for  MACE  behaved 
very  closely  to  zero  percent  damping  as  the  results  suggest. 


5.4.5  Why  shouldn’t  it  work  so  well? 

We  have  reasoned  why  the  experimental  data  should  follow  the 
theoretical  data,  but  we  can  also  reflect  on  why  it  shouldn’t.  Although  the 
plant  is  not  moving  very  much,  it  is  still  a  real-life  nonlinear  system  with 
all  its  unpredictable  characteristics.  There  are  many  real  hardware 
hurdles  that  this  test  overcomes:  there  is  fiiction  and  stiction  in  the 
actuator  and  in  the  joints  of  the  bus;  the  computer,  amplifiers,  and  filters 
all  have  dynamics  of  their  own;  energy  migrates  between  frequencies  as  in 
all  real  systems;  and  the  MACE  suspension  system  adds  funny  dynamics 
to  the  structure  including  the  plimge  modes,  the  pendulum  modes,  and  the 
violin  modes  of  the  rods.  All  of  these  factors  are  not  accounted  for  in  the 
theoretical  formulation.  Although  the  test  experiences  all  of  these  “non- 
linearities,”  the  frequency  response  data  still  follows  closely  what  we  would 
predict  firom  a  simple  linear  model.  We  will  see  in  the  next  section  that 
once  we  start  slewing  large  angles,  this  linearity  assumption  becomes  less 
valid. 


5.5  Closed  Loop  Time  Histories 

5.5.1  Slew  profile  and  limitations 

All  the  slew  calculations  were  done  in  real  time  with  user  inputs.  The 
slew  profile  was  originally  derived  firom  scaled  slews  of  what  a  larger  size 
pointing  payload  would  do,  such  as  a  telescope  on  the  EOS  platform.  This 
was  determined  to  scale  to  an  angle  of  120  deg  (±60**),  a  velocity  of  50  d^sec, 
and  an  acceleration  of  125  de^s^.  The  slew  angle  of  the  gimbal  is  governed 
in  actual  hardware,  however,  by  the  limitations  of  the  suspension  system. 
The  suspension  system  travel  does  not  allow  for  any  slews  larger  than  20 
degrees  as  the  pressure  pistons  bottom  out  for  larger  moves.  The  move  was 
then  transformed  for  a  20  deg  slew  that  kept  that  same  maximum  velocity. 
The  slew  has  a  trapezoidal  velocity  profile  which  has  a  1-2>1  ratio  of 
acceleration  -  constant  velocity  •  deceleration.  The  profile  in  shown  in 
figure  5.10.  This  profile  was  used  as  a  typical  slew  for  a  pointing  payload  in 
our  experiments.  The  slew  formulation  was  programmed  into  the  real¬ 
time  computer  to  convert  a  user  step  input  command  into  this  type  of 
profile.  Since  this  smoothing  function  manipulates  inputs,  it  can  also  be 
considered  a  type  of  real-time  input  shaping.  It  does,  though  not  very 
efficiently,  filter  out  unwanted  high  firequency  content  firom  the  input. 

5^5.2  Unshaped  response 

A  siunmaiy  of  the  unshaped  20  degree  slew  is  shown  in  figures  5.11  - 
5.13.  We  can  see  that  the  controller  follows  the  move  fidrly  closely  with  a  bit 
of  time  lag  and  some  steady  state  error,  as  we  would  expect  from  a  PD 
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SLEW  PROFILE 


FIGURE  5.10:  Profile  of  gimbal  an^e  command 


GIMBAL  ANGLE  RESPONSE  TO  20“  UNSHAPED  SLEW 
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FIGURE5.il:  Gimbal  angle  command  and  response  for  a 
smooth  20**  slew. 
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controller.  The  slew  is  smooth  enough  such  that  there  is  no  visible 
overshoot  or  settling  time  in  the  gimbal  angle  response.  The  velocity  plot 
shows  the  trapezoidal  profile  that  the  gimbal  is  following.  There  is  a  bit  of 
overshoot  in  the  velocity  curve  which  shows  up  as  a  dip  at  the  top  of  the 
trapezoid.  The  maximum  velocity  is  around  50  deg/sec.  The  accelerometer 
output  shows  the  response  at  the  other  end  of  the  bus.  The  acceleration  plot 
shows  the  multiple  mode  behavior  of  the  response,  and  that  most  of  the 
amplitude  is  governed  by  the  first  bending  mode  of  the  system  (1.75Hz).  We 
can  also  note  from  the  plot  the  slight  0.2Hz  plimge  mode  of  the  susi)ension 
system.  This  will  become  more  evident  in  the  shaped  responses.  As 
discussed  in  Chapter  3,  this  mode  is  an  order  of  magnitude  below  the  first 
resonance  of  the  bus  and  it  should  have  little  effect  on  the  dynamics  of  the 
plant.  The  vibration  in  the  acceleration  response  is  large  and  lasts  for  quite 
a  while,  begging  the  need  for  a  vibration  control  sdtieme. 


5.5.3  Shaped  Response 

5.5.3.1  One  mode  shaped 

Figures  5.14  -  5.15  show  a  summary  of  what  the  slew  looks  like  when  it  is 
convolved  with  a  three  impulse  shaper  for  the  first  bending  mode.  A  close 
look  at  the  slew  profile  reveals  that  the  shaped  slew  reaches  the  end  position 
about  .57  seconds  (l/1.75Hz)  later  than  the  unshaped  move  in  figure  5.11. 
This  is  the  time  delay  cost  or  the  length  of  the  shaper.  The  acceleration 
response  exhibits  a  significant  reduction  in  the  vibration  amplitude.  The 
first  mode,  which  is  the  most  visible  and  the  longest  lasting,  is  almost 
unnoticeable.  In  the  nine  seconds  of  data  that  is  taken,  we  can  see  most  of 
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fRE5.14:  Gimbal  angle  command  and  response  for  a  20® 
slew  shaped  for  one  mode  (1.75Hz). 


ACCELERATION  RESPONSE  TO  SHAPED  20“  SLEW  (1  MODE) 
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Node  4  vertical  acceleration  response  to  20®  slew 
shaped  for  one  mode  (l.TSHz). 
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FAST  FOURIER  IRANSFORM  OF  SHAPED  20”  SLEW  (1  MODE) 


FRE9JENCY(Hz) 


FIGURE  5.16;  Fast  Fourier  Transfom  of  residual  vibration  at  node 
4  acceleroDseter  for  20®  slew  shaped  for  Imode  (1.75Hz). 

the  vibration  in  the  other  modes  settle  as  well.  A  more  quantitative  way  to 
view  this  data  is  to  take  an  FFT  (Fast  Fourier  Transform)  of  the  data  to  get  a 
frequency  spectrum  over  a  range  of  frequencies.  To  measure  the  residual 
vibration,  we  need  only  to  FFT  the  data  after  the  move  has  been  completed, 
which  in  this  case  is  at  about  l.lsec  (.53-I-.57).  Five  runs  were  taken  and 
averaged  to  plot  the  data  shown  in  figure  5.16.  The  amplitude  of  the 
vibration  is  plotted  on  a  linear-linear  scale  for  the  0-20Hz  range.  We  can  see 
the  reduction  of  the  peak  at  1.75Hz.  This  corresponds  to  an  effective 
vibration  reduction  of  97%.  We  also  see  that  substantial  vibration  has  been 
canceled  in  the  third,  fifth,  and  sixth  modes.  None  of  those  modes  were 
canceled  to  a  satisfactory  level,  but  it  is  interesting  to  note  the  consequence 
of  the  repeating  nature  of  the  impulse  shaper.  This  shaper  is  a  very 
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effective  one  for  the  first  mode  and  correlates  closely  to  the  results  shown  in 
the  frequency  response  experiments.  We  see  in  the  next  section  that  the 
higher  modes  are  more  difficult  to  suppress. 


5.S.3.2  Four  modes  shaped 

We  now  move  on  to  the  data  taken  with  the  four  mode  direct  solution 
impulse  shaper.  The  same  plots  as  in  the  previous  section  are  shown  in 
figures  5.17  -  5.18.  We  notice  that  the  cost  in  shaper  delay  is  slightly  higher 
(around  .676  seconds).  The  acceleration  response  plot  shows  an 
improvement  over  the  one  mode  sequence.  This  is  a  slower  move  than  the 
one  mode  sequence,  so  we  see  less  vibration  during  the  move  also.  Almost 
all  of  the  visible  vibration  is  gone  by  8  seconds  into  the  move.  The  FFT  data 
in  figure  5.19  shows  in  more  detail  the  improvement  over  the  unshaped 
move  and  the  one  mode  shaped  move.  The  reduction  in  the  first  mode  is 
practically  unchanged.  The  vibration  reduction  is  about  96%.  The  peaks  at 
8.9Hz  and  9.13Hz  have  been  further  reduced.  However,  their  amplitudes  are 
still  very  visible.  In  the  time  traces,  we  can  see  that  the  move  has  ended  at 
about  1.2  seconds,  but  there  is  still  very  visible  vibration  with  a  frequency  of 
around  9Hz.  The  percentage  reductions  here  were  only  69%  and  85% 
respectively.  This  is  inconsistent  with  the  theoretical  and  experimental 
insensitivity  curves  which  showed  almost  zero  percent  vibration  over  the 
whole  range  between  8  and  lOHz.  These  are  highly  excitable  modes  and 
that  may  contribute  to  their  sensitivity,  but  we  would  expect  the  robustness 
of  this  4  mode  sequence  to  overcome  any  tendency  for  those  modes  to  vibrate. 
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In  the  next  section,  we  try  using  more  robust  shapers  to  cancel  the  vibration 
in  this  region. 


5.6  Using  More  Robust  Sequences 

We  need  now  to  address  why  the  4  mode  direct  solution  shaper  is  not  as 
effective  for  the  higher  modes  at  8.9  and  9.13Hz  as  it  was  for  the  1.75Hz 
mode.  The  problem  lies  in  the  fact  that  impulse  shaping  is  based  on  linear 
theory  and  here  it  is  being  applied  to  a  nonlinear  system.  We  saw  hints  of 
linearity  in  the  frequency  response  data  and  we  may  expect  that  the  time 
histories  would  be  somewhat  dose  to  that,  if  not  the  same.  But  it  is  not  The 
movement  of  the  gimbal  over  a  large  angle  induces  unmodeled  effects  that 


GIMBAL  ANGLE  RESPONSE  TO  20°  SHAPED  SLEW  (4  MOIXES) 


FIGURE  5.17:  Gimbal  angle  command  and  response  to  20**  slew  shtqjed 
for  4  modes  (1.75,  5.23,  8.68,  9.13Hz). 
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disable  tbe  impulse  shaper  from  canceling  higher  frequency  modes.  We 
will  from  here  on  call  these  unwanted  effects  “nonlinearities.”  We  now 
want  to  determine  if  a  more  robust  sequence  can  overcome  these  non- 
linearities  and  cancel  the  remaining  residual  vibration  in  the  8.9Hz  and 
9.13Hz  modes. 

As  mentioned  in  Chapter  2,  the  direct  solution  method  loses  some 
robustness  in  higher  modes  when  compared  to  the  convolution  method. 
Seeking  more  robustness  in  the  8  to  lOHz  region,  a  convolved  sequence  of  the 
two  modes  around  9Hz  was  formulated.  Four  impulse  sequences  were  also 
considered,  adding  another  derivative  constraint  to  the  zero  vibration 
eqxiations.  Two  new  four  impulse  sequences  were  calculated;  one  sequence 
for  the  frequency  in  between  the  two  modes  (8.9Hz),  and  the  second,  a 
convolved  sequence  for  both  modes.  Both  convolution  and  second  derivative 
constraints  add  substantially  to  the  length  of  the  sequences.  For  example,  a 
direct  solution  sequence  using  two  impulses  per  mode  for  the  two 
frequencies  at  8.68  and  9.13Hz  yielded  a  sequence  0.112  seconds  long  where  a 
convolved  four  impulse  sequence  is  0.337  seconds  long.  What  we  lose  in 
time  delay,  we  are  gaining  in  robustness.  The  new  sequences  concentrate 
on  the  modes  in  the  9Hz  area,  foregoing  canceling  vibration  at  the  other 
modes.  An  exhaustive  list  of  impiilse  sequences  were  implemented  in  an 
attempt  to  cancel  the  vibration  at  these  frequencies.  However,  none  of  them 
were  successful.  The  results  of  the  slews  with  the  different  shapers  are 
shown  in  figures  5.20  -  5.21.  All  of  these  plots  are  averages  of  five  runs  and 
contain  164  points  in  the  0-20Hz  region.  This  plots  are  not  intended  to  show 
the  exact  amplitudes  of  the  individual  shaped  moves.  There  is  too  much 
data  here  to  show.  It  does  reveal  the  trends  in  the  data;  that  no  matter  what 


shaper  was  used,  the  response  was  always  pretty  much  the  same  in  that 
region.  Even  the  convolved  four  impulse  sequence,  which  contains  16 
impulses  and  is  extremely  robust  in  the  8-lOHz  range,  foiled  to  effectively 
cancel  those  modes.  Interestingly,  the  shai)er  we  would  least  expect  to  be 
effective,  the  three  impulse  sequence  for  the  first  mode,  performed  about 
median  among  all  the  other  shapers. 

The  results  show  that  a  more  robust  sequence  cannot  cancel  the 
remaining  30%  or  so  vibration  left  in  the  fifth  and  sixth  modes.  In  fact, 
there  was  very  little  deviation  in  the  residual  vibration  for  all  the  shapers 
that  were  implemented  It  follows  to  reason  that  the  resulting  vibration  is 
caused  by  something  other  than  what  impulse  shaping  theory  would 
predict.  This  “nonlinear”  effect  is  evident  in  the  higher  fii^quendes  when 
we  slew  large  angles.  From  the  results  of  these  experiments,  we  may  say 
that  nothing  in  the  domain  of  impulse  shaping  seems  to  be  able  to  overcome 
the  nonlinearities  associated  with  large  angle  movement  to  effectively 
suppress  vibration  in  the  higher  fii^uency  modes. 
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It  is  difficult  to  determine  which  particular  nonlinearity  causes  most  of 
the  vibration  in  this  region.  Most  likely,  it  is  a  combination  of  many 
different  things.  We  have  identified  some  nonlinearities  as  being  small 
from  the  success  of  the  frequency  response  experiments.  For  example; 
amplifier/computer/filter  dynamics,  friction,  stiction,  digitization,  and  the 
suspension  effects  were  all  found  to  be  negligible  factors.  Of  course,  the 
slew  experiments  induce  different  amplitudes  of  these  factors,  but  we 
wouldn’t  expect  any  one  of  them  to  be  the  major  contributing  factor  to  the 
residual  vibration. 

Some  other  nonlinearities  can  be  examined  as  well.  Centripetal 
acceleration  is  caused  when  the  gimbal  slews  large  angles.  This,  however, 
is  found  to  be  very  small.  For  a  maximum  velocity  of  50  deg/sec,  this 
translates  to  a  centrifugal  force  of  less  then  .007  Ibf.  And  this  is  applied  for 
only  a  fraction  of  the  slew  time.  The  effect  on  the  dynamics  of  the  system  is 
minimal.  There  is  also  a  kinematic  nonlinearity  of  frequencies  shifting  as 
the  configuration  of  the  plant  changes  as  discussed  in  Chapter  4.  In 
hardware,  this  was  found  to  be  small,  as  a  transfer  function  of  the  plant 
with  the  gimbal  at  20°  showed  less  than  1%  change  in  all  the  firequendes 
under  lOHz.  The  gimbal  is  small  in  mass  and  in  length  compared  to  the 
bus,  so  its  movement  causes  little  change  in  the  firequency  characteristics  of 
the  whole  structure.  The  gravity  torque  was  also  an  issue.  This  is  a  very 
important  issue  in  terms  of  space-home  systems  such  as  MACE,  since 
gravity  is  a  factor  only  during  the  ground  testing  phase.  The  gravity  is 
pretty  linftgr  in  the  ±20°  range  however.  Its  effect  on  the  motor  was  found  to 
be  linear  to  ±2%  in  this  region.  Moreover,  the  indusion  of  a  gravity 
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compensator,  as  mentioned  in  Chapter  3,  did  not  remedy  the  problem.  All 
of  these  factors  were  considered  as  contributors  to  the  remaining  vibration, 
but  not  one  of  them  was  found  to  be  nuyor. 

The  severity  of  the  nonlinearity  does  seem  to  be  a  function  of  the  size  of 
the  slew.  Different  size  slews  of  the  same  maximum  velocity  were  used  in  a 
series  of  experimental  runs  with  the  same  impulse  shaper.  The  slews  have 
been  scaled  in  both  time  and  angle  to  retain  the  same  velocity  profile  for  the 
smaller  an^e  dews.  So,  a  slew  to  10**  would  have  the  same  shape  as  the  20** 
slew  but  would  finish  its  move  in  half  the  time.  A  plot  of  the  amplitude  of 
vibration  in  the  8-lOHz  range  for  the  various  angle  slews  with  the  two  mode 
three  impulse  convolved  sequence  is  shown  in  figure  5.22.  We  can  see  that 
the  residual  vibration  is  more  severe  with  the  20^*  move  and  is  the  least 
visible  in  the  move.  This  is  despite  the  fact  that  the  unshaped  smaller 
angle  moves  are  actually  putting  more  energy  into  these  modes  because 
they  are  more  abrupt,  sharper  moves.  This  was  done  to  accent  the  higher 
frequency  modes  for  the  smaller  angle  slews.  Small  angle  slews  done  in 
the  same  time  as  the  20°  run  would  have  caused  little  vibration  and 
provided  little  insight  into  the  nonlinearity  of  large  angle  slews. 

The  data  suggests  that  the  nonlinear  effect  is  linked  to  the  size  of  the 
move.  As  mentioned  in  the  frequency  response  analysis  in  section  5.4, 
linear  approximations  are  generally  valid  in  systems  with  small 
amplitudes.  Larger  amplitudes  tend  to  cause  to  system  to  act  more 
nonlinearly  and  thus  less  responsive  to  impulse  shaping.  The  slews  done 
here  have  much  higher  amplitudes  than  the  frequency  response 
experiments.  More  energy  from  the  amplifiers  causes  the  structure  to 
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L23:  Fast  Fourier  Transfonn  of  residual  vibration  amplitude 
for  1**  slew  shaped  for  2  modes  (8.68, 9.13Hz). 
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physically  deform  more.  What  we  are  observing  is  a  system  which  acts 
linearly  only  for  small  angle  displacements  (<!**). 

An  example  of  a  small  angle  move  is  shown  in  figure  5.23.  This  1°  move 
shows  similar  effectiveness  as  that  was  seen  in  the  frequency  response 
plots.  The  move  is  short  and  abrupt  and  causes  a  lot  of  amplitude  in  the  8- 
lOHz  range.  The  shaped  response  cuts  that  amplitude  by  99%.  The 
resulting  vibration  is  even  less  than  that  left  by  the  20**  slew  which  had 
much  less  unshaped  vibration  amplitude.  Although  the  20**  slew  puts  less 
energy  into  the  higher  modes,  it  puts  more  total  energy  into  the  system  than 
the  smaller  slews  and  causes  much  more  amplitude  during  the  move  time. 
The  larger  amplitude  of  vibration  and  the  greater  change  in  plant 
configuration  contribute  to  the  nonlinearity  of  the  system  and  cause  the 
impulse  shaper  to  be  less  effective. 

We  can  conclude  that  the  a  good  portion  of  the  nonlinear  effect  is  caused 
by  the  large  vibration  amplitudes  and  change  in  configuration  due  to  large 
angle  slews.  Many  other  nonlinearities  were  explored  as  possibilities,  but 
not  one  was  found  to  be  a  m^or  factor.  It  is  likely  that  a  combination  of  all 
of  them  contribute  to  the  nonlinearity  of  the  whole  system.  The  size  of  the 
move  was  seen  to  have  the  most  correlation  to  the  amount  of  vibration  that 
was  canceled  by  the  impulse  shaper. 
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5.8.1  Breaking  with  impulse  shaping  assumptions 

Impulse  shaping  only  assumes  the  plant  eigenvalues  in  the  system  it  is 
trying  to  control.  The  plant  is  seen  as  a  black  box  with  some  unique  second 
order  poles  associated  with  it.  The  type  of  controller  determines  where 
these  poles  are,  and  in  some  cases  can  even  add  poles  to  the  system.  An 
impulse  shaper  is  concerned  only  with  these  closed  loop  poles  and  tries  to 
minimize  the  energy  put  into  those  frequencies.  The  black  box,  however, 
has  components  which  cannot  be  ignored.  In  this  section,  we  will  show 
that  the  controller  inside  the  system  is  an  important  factor  in  determining 
the  efiTectiveness  of  an  impulse  shaper. 

A  block  diagram  of  MACE  may  look  like  the  following: 


SYSTEM 


FIGURE  5.24:  A  system  diagram  of  MACE 


The  system  that  the  shaper  is  tr3dng  to  affect  consists  of  the  gimbal 
actuator,  its  controller,  and  the  bus.  The  gimbal  and  bus  can  be  seen  as 
decoupled  here  because  of  the  large  relative  size  of  the  bus  and  the 
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configuration  of  the  gimbal  sensor  (the  encoder).  The  gimbal  as  the 
actuator  has  much  authority  to  induce  vibration  in  the  bus,  but  there  is  not 
much  feedback  from  the  bus  to  the  gimbal  angle.  Of  course,  both 
components  are  part  of  the  same  system  and  the  same  resonances  are  seen 
on  both  components,  but  the  large  vibration  that  the  gimbal  can  cause  in  the 
bus  is  not  significant  in  the  angle  amplitude  of  the  gimbal.  This  is  due  to 
the  fact  that  the  encoder  is  an  angle  sensor  where  most  of  the  vibration 
induced  in  the  bus  is  vertical.  This  is  also  partly  the  reason  a  PD  loop  was 
able  to  be  successfully  closed  on  the  gimbal. 

The  system  configuration  is  one  closed  loop  system  (gimbal)  which 
excites  an  open  loop  system  (bxis).  The  sensor  that  we’re  trying  to  control 
(the  accelerometer)  is  not  a  part  of  the  closed  loop  system.  This  presents  an 
interesting  challenge  in  that  the  controller  now  becomes  an  integral  part  of 
how  much  vibration  is  induced  in  the  bus.  The  torque  that  the  controller 
applies  at  the  gimbal  is  in  effect  the  input  into  the  plant  since  the  majority 
of  the  dynamics  is  in  the  bus.  In  impulse  shaping  theory,  this  can  still  be 
seen  as  a  black  box  with  unique  resonances.  The  controller  only  determines 
where  and  how  large  the  poles  are.  However,  the  controller  is  important  in 
terms  of  impiilse  shaping  because  it  determines  how  much  authority  the 
actuator  has  over  the  dynamics  of  the  bus. 

Controllers  tend  to  act  like  low-pass  filters  in  that  the  gain  rolls  off  past 
the  bandwidth  frequency  of  the  controller  so  that  the  high  frequency 
amplitudes  are  minimal.  For  example,  the  current  PD  controller  of  the 
gimbal  could  not  follow  a  20Hz  sine  wave  very  well.  It  woiild  have  trouble 
with  any  move  that  has  frequency  content  greater  than  its  bandwidth.  This 
is  the  reason  that  the  controller  is  usually  designed  with  the  types  of  likely 
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command  inputs  in  mind.  Similarly,  a  low  bandwidth  controller  will  not  be 
able  to  follow  the  precise  and  quickly  changing  path  demanded  by  a 
complicated  impulse  shaper.  This  will  be  more  true  when  higher  modes 
are  being  shaped  and  the  impulses  are  closer  together.  For  the  shaper  to  be 
most  effective,  the  path  of  the  command  needs  to  be  followed  closely.  This 
ensures  that  the  vibration  induced  during  the  move  will  be  canceled  by  the 
impulses  of  the  shaper.  Also  the  move  needs  to  be  consistent.  Deviations  at 
irregular  times  during  the  move  can  render  the  shaper  ineffective.  This  is 
the  reason  the  bandwidth  of  a  controller  is  a  factor  in  how  effective  a  shaper 
can  be.  One  may  argue  that  a  low  bandwidth  controller  won’t  put  high 
frequency  energy  into  the  system  in  the  first  place,  so  that  canceling 
vibration  at  those  frequencies  is  not  an  issue.  If  the  gain  at  that  frequency 
were  zero,  then  this  would  be  true.  However,  this  is  never  the  case,  and 
modes  that  have  veiy  high  gains  do  not  need  much  excitation  to  cause  a  lot 
of  vibration.  We  will  see  in  the  next  section  that  for  highly  excitable  modes, 
such  as  the  fifth  and  sixth  modes  of  MACE,  the  vibration  is  still  an  issue 
even  for  low  bandwidth  controllers. 

The  bandwidth  of  the  controller  determines  how  closely  an  actuator  can 
follow  a  command  and  determines  the  consistency  of  the  move.  A  higher 
bandwidth  controller  will  be  more  able  to  follow  high  firequency  moves  since 
its  frequency  response  will  be  higher  in  those  frequency  ranges.  Also  it  will 
be  more  able  to  overcome  inconsistencies  of  the  system  (friction, 
disturbances,  etc.).  This  makes  for  a  more  robtist  and  consistent  move. 
This  is  also  very  important  in  the  precise  implementation  of  impulse 
sequences.  It  follows,  then,  that  the  bandwidth  of  the  gimbal  controller  can 
determine  how  effective  a  shaper  is  in  reducing  residual  vibration.  Since 
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the  gimbal  is  the  cause  of  the  vibration  in  the  bus,  a  stifier  gimbal  controller 
will  more  precisely  apply  inputs  to  the  bus.  Because  the  impulses  require 
very  precise  amplitudes,  this  determines  then  how  much  vibration  we  can 
cancel  in  the  higher  modes.  We  will  now  look  at  results  firom  shaping 
experiments  done  with  different  controllers. 


5.8.2  Designing  different  bandwidth  controllers 

The  controller  used  in  the  experiments  thus  far  was  designed  to  have  a 
bandwidth  of  lOHz  as  mounted  to  a  rigid  frame.  This  means  that  the  gain 
of  the  closed  loop  transfer  function  is  .707  at  lOHz.  Bandwidth  becomes  an 
ill-defined  term  for  a  PD  controller  once  the  gimbal  is  attached  to  the  rest  of 
MACE  because  there  are  now  unmodeled  modes  which  are  within  the 
bandwidth.  The  controller  was  designed  for  the  gimbal  itself,  not  for  the 
gimbal  mounted  to  the  bus  assembly.  If  we  look  on  the  transfer  function 
plot  for  the  lOHz  bandwidth  controller  (figure  5.25),  we  see  the  peaks  and 
valleys  before  the  lOHz  point  In  this  case,  the  closed  loop  transfer  function 
crosses  the  .707  point  more  than  once,  so  it  is  difficult  to  define  bandwidth 
here.  Moreover,  the  dynamics  of  the  system  have  changed  such  that  the 
previous  .707  crossing  point  at  lOHz  has  shifted  due  to  the  addition  of  a  zero. 
Thus,  we  will  use  the  term  bandwidth  as  a  measurement  of  the  relative 
strength  or  stiffiiess  of  the  controllers. 

Both  a  3Hz  bandwidth  and  20Hz  bandwidth  controller  were  designed  in 
addition  to  the  old  lOHz  bandwidth  controller.  A  summary  of  the 
parameters  of  the  3  existing  controllers  is  shown  in  table  5.1. 
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All  three  are  basically  PD  controllers  with  some  stabilizing  features  in 
the  loop.  The  3Hz  bandwidth  controller  did  not  require  a  delay  buffer  as 
there  was  enough  gain  and  phase  margin  to  account  for  the  delay  of  the 
computer  processor  board.  This  controller  is  very  weak  but  also  very  stable. 
The  20Hz  controller  employed  the  delay  buffer  and  also  required  a  notch 
filter  in  the  loop.  A  mode  at  85Hz  crossed  the  gain  margin  in  closing  the 
loop  and  went  unstable.  A  second  order  Butterworth  notch  filter  with 
comer  frequencies  at  80Hz  and  90Hz  proved  to  be  an  effective  solution  to  the 
problem.  This  controller  is  stable  and  very  stiff.  The  closed  loop  transfer 
functions  from  the  angle  command  to  the  encoder  angle  output  are  shown 
in  figure  5.25  for  all  three  controllers. 
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CLOSED  LOOP  TRANSFER  FUNCTIONS  FOR  3  CONTROLLERS 


FIGURE  5.25;  Cosed  loop  transfer  functions  of  ccnmnanded  angle 
to  encoder  output  for  three  different  ccmtrollers. 

5.8.3  Effectiveness  of  shaping  on  different  controllers 

As  we  have  discussed,  the  bandwidth  of  the  controller  can  determine 
how  effective  a  shaper  is  in  canceling  vibration.  This  is  because  the  higher 
bandwidth  controller  follows  the  move  more  closely  and  is  more  true  to  the 
command  that  the  shaper  demands.  The  same  inputs,  unshaped  and 
shaped,  were  employed  with  all  three  controllers.  The  original  .53  second 
20^’  slew  was  used.  Five  runs  of  the  data  were  recorded  for  each  controller, 
and  the  Fast  Fourier  Transform  of  the  residual  vibration  was  taken  for 
each  and  averaged.  The  results  for  the  three  controllers  are  shown  in 
figures  5.26  -  5.28.  The  lOHz  controller  response  is  one  weVe  already  seen. 
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[GURE  5.27:  Qoseup  view  of  FFT  of  residual  vibration  for  3Hz 
controller.  20**  stew  is  shaped  f(V  two  modes  (8.68, 9.13Hz). 
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FIGURE  5.28:  Qosetq)  view  <tfFFTdr  residual  vibratkxi  for  20Hz 
coatxoller.  20**  slew  is  shaped  for  two  inodes  (8.68, 9.13Hz). 


COMPARISON  OF  DIFFERENT  BANDWIDTH  CONTROLLERS 


FREQUENCY  (Hz) 

FIGURE  5.29:  A  cotnpariscn  of  residual  vibiati(maoq;>litudes  in 
die  6  to  12Hz  region  for  three  different  (x»tndlers.  Same  20**  slew 
is  shaped  for  two  modes  in  each  case  (8.68, 9.13Hz) 
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It  is  interesting  to  note  that,  in  using  the  SELz  bandwidth  controller,  we 
actually  incur  more  vibration  in  the  8-lOHz  range  when  the  move  is  shaped. 
This  is  a  case  where  the  controller,  as  seen,  has  very  little  authority  to  put 
energy  into  these  modes  whether  shaped  or  unshaped.  The  frequency 
response  of  the  controlled  plant  is  very  low  in  this  range.  However,  the 
amplitude  of  these  modes  are  very  high,  and  are  nevertheless  ezdted  by  the 
3Hz  controller.  The  controller  has  a  hard  time  trying  to  follow  the  shaped 
move  because  of  its  low  feedback  gains.  In  fact,  for  a  20**  command,  the 
controller  can  only  push  the  gimbal  up  to  a  steady  state  angle  of  16**. 
Therefore,  the  unshaped  move  shows  little  vibration  in  this  range.  It  is  a 
noteworthy  fact,  however,  that  the  shaped  move  has  a  detrimental  effect  on 
the  vibration  cancellation.  We  Min  conclude  that  the  controller  could  not 
follow  the  move  dosely  enough  to  have  any  impact  on  those  modes  at  all. 
The  extra  vibration  comes  from  the  fact  that  the  move  is  longer  and  more 
inconsistent.  This  gives  more  chance  for  system  disturbances  or 
inconsistendes  to  add  hifd^  firequency  energy  to  the  resulting  output. 

The  20Hz  bandwidth  controller  shows  more  promising  results.  As 
shown  in  the  figure,  the  fifth  and  sixth  modes  have  higher  imshaped 
amplitudes  than  the  other  controllers.  This  is  due  to  the  same  reason  that  a 
shaped  move  puts  less  energy  into  higher  modes.  The  move  is  quicker, 
crisper,  and  faster  such  that  the  profile  has  more  high  frequency 
components.  The  firequency  response  of  the  controller  in  this  range  is 
nearly  1  for  this  controller,  so  it  can  follow  the  shaped  moves  fairly  dosely. 
The  shaped  amplitudes  are  about  7%  and  14%  of  the  unshaped  amplitudes 
in  the  fifth  and  sixth  modes.  This  is  still  not  the  level  of  performance  that 
we  saw  on  the  frequency  response  experiments,  but  it  constitutes  a 
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significant  improvement  over  the  lower  bandwidth  controllers.  Although 
the  controller  is  more  true  to  the  demanded  shaped  move,  there  are  stiU  the 
nonlinearities  associated  with  a  20^  slew  such  that  all  the  vibration  cannot 
be  canceled. 

Not  only  is  the  higher  bandwidth  controller  an  improvement  in  terms  of 
percentage  vibration  reduction,  but  it  also  has  the  absolute  least  residual 
vibration.  A  close-up  of  the  8-lOHz  region  for  the  FFTs  of  the  shaped 
responses  of  the  three  controllers  is  shown  in  figure  5.29.  The  plot  reveals 
that,  even  in  absolute  terms,  the  20Hz  controller  performed  the  best  out  of 
the  three.  This  is  despite  the  fact  that  the  unshaped  slew  caused  much 
more  vibration.  The  20Hz  controller  is  a  win  in  both  percentage 
cancellation  and  absolute  residual  amplitude. 


5.8.4  Conclusions 

We  have  just  touched  on  a  topic  that  has  been  mostly  ignored  in  past 
research  on  input  shapers.  That  is  the  role  of  the  controller  in  determining 
the  effectiveness  of  a  shaping  algorithm.  We  have  seen  here  through 
experimental  results  that  the  controller  does  determine  how  effective  an 
impulse  shaper  can  be  in  this  application.  The  20Hz  bandwidth  controller 
was  shown  to  be  more  effective  in  canceling  vibration  in  the  8-lOHz  range 
than  either  the  3Hz  or  lOHz  bandwidth  controllers.  The  3Hz  controller  even 
showed  a  slight  detrimental  effect  on  the  shaped  response.  Thus  far  we 
have  only  dealt  with  bandwidth  of  PD  controllers.  Most  likely,  higher  order 
controUers  will  also  have  an  effect  on  the  shaper  effectiveness.  The  more 
consistently  and  precisely  a  shaped  command  is  input  into  the  system,  the 


more  efitectiye  an  impulse  shaper  can  be.  A  more  sophisticated  pointing 
controller  for  the  MACE  gimbal  should  show  an  improvement  in  vibration 
cancellation. 

5.9  Closure 

In  this  chapter,  we  have  looked  at  various  hardware  experiments  to 
determine  the  effect  of  an  impulse  shaper  on  the  residual  vibration  of  the 
MACE  test  article.  We  began  by  defining  the  performance  metric.  We  then 
went  on  to  look  at  firequency  response  data  with  the  shaper  in  the  system. 
These  proved  to  be  very  effective  in  canceling  unwanted  vibration  and 
followed  the  theoretical  data  very  closely.  The  time  histories  of  slews  to  20*^ 
showed  less  success  in  percentage  vibration  reduction.  The  first  bending 
mode  was  effectively  canceled,  but  the  higher  modes,  specifically  the  fifth 
and  sixth  modes  between  8  and  lOHz  showed  significant  amplitude  even  in 
the  shaped  runs.  A  look  at  more  robust  impulse  shapers  showed  that 
nothing  in  the  realm  of  linear  impulse  shaping  seemed  to  be  able  to 
overcome  the  nonlinearity  associated  with  a  large  angle  slew.  The  most 
simple  impulse  shaper  was  just  as  effective  as  any  of  the  more  complicated 
and  more  robust  sequences.  Slews  done  for  different  angles  at  the  same 
TtiflTiTniim  velocity  showed  that  this  nonlinearity  is  a  function  of  the  size  of 
the  slew.  For  a  one  degree  slew,  the  response  looked  similar  to  that  of  the 
transfer  function  data. 

We  then  looked  at  how  the  controller  in  the  S3rstem  can  change  the 
effectiveness  of  an  impulse  shaper.  Three  different  bandwidth  controllers 
were  designed  and  tested  with  the  same  set  of  slews.  The  results  showed 
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that  a  higher  bandwidth  controller  is  more  effective  in  fianpAling  vibration 
at  higher  frequencies.  This  is  because  higher  bandwidth  controllers  will 
more  accurately  and  consistently  follow  the  move  demanded  by  the  impulse 
shaper,  which  tends  to  be  more  irregular  and  have  high  frequency  content. 
The  20Hz  bandwidth  controller  was  found  to  be  most  effective  in  percentage 
reduction  and  even  in  absolute  terms  in  leaving  the  least  amount  of 
residual  vibration. 


Chapter  6 


6.1  Conclusions 

We  have  investigated  the  issues  involved  with  implementing  impulse 
command  shapers  on  a  large  flexible  space  structure.  There  was  a  need  to  do 
this  because  of  a  lack  of  real  experimental  data  to  verify  the  effectiveness  of 
impulse  shaping  on  complex  multiple  mode  systems.  My  contribution  to  the 
field  of  input  shaping  has  been  to  answer  some  questions  about  which 
qualities  of  impulse  shaping  work  well  and  which  ones  don’t.  I  also  suggest 
ways  in  which  to  make  shaping  more  effective  in  implementation. 

In  Chapter  2,  we  discussed  the  theoretical  backgroimd  of  the  impulse 
sequences.  The  zero  vibration  constraint  equations  are  derived  firom  the 
response  of  second  order  systems  to  impulses.  Derivative  constraints  are 
added  to  give  the  sequence  more  robustness  to  variations  in  frequency  and 
damping  ratio.  We  discussed  methods  of  generating  sequences  for  multiple 
mode  systems  and  contemplated  the  tradeoffs  between  the  convolution 
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method  and  the  direct  solution  method.  Some  suggestinas  were  then  made  on 
finding  the  most  efficient  sequence  for  particular  applications.  The  variables 
that  were  considered  are:  the  number  of  modes  to  shape,  direct  solution 
versus  convolution,  and  the  number  of  impulses  per  mode.  These  can  be  used 
in  our  favor  to  find  the  sequence  that  best  suits  a  particular  application. 

The  testbed  for  the  experiments,  MACE,  was  discussed  in  Chapter  3. 
MACE  is  a  NASA/MIT  active  control  experiment  that  is  scheduled  to  fly  on 
the  space  shuttle  in  the  summer  of  1994.  Its  goal  is  to  provide  a  platform  for 
investigating  zero  gravity  dynamics  of  a  large  space  structure  which  utilizes 
Controlled  Structures  Technology.  A  description  of  the  Development  Model 
hardware  was  provided.  Some  results  firom  d}mamic  testing  were  presented 
and  a  modal  description  of  the  test  article  was  presented.  We  then  went  on  to 
detail  some  of  the  controller  issues  in  the  gimbal  actuator  ar  1  the  AC- 100 
computer  system.  The  MACE  test  article  was  found  to  be  a  very  suitable 
testbed  for  the  goals  of  this  research. 

Non-linear  simulation  was  the  topic  of  Chapter  4.  DISCOS  was  used  as  a 
tool  for  verif^dng  the  effectiveness  of  input  shaping  on  both  open  loop  and 
closed  loop  configurations  of  MACE.  The  kinematic  non-linearity  due  to 
changes  in  the  configuration  of  the  plant  was  discussed.  The  frequency 
change  in  the  structure  was  shown  to  be  a  hindrance  in  the  canceling  of  the 
first  mode. 

The  experimental  data  ^as  presented  in  Chapter  5.  It  was  found  that  in 
frequency  response  tests,  impulse  shaping  worked  nearly  to  the  level 
predicted  by  theory.  Insensitivity  curves  for  both  experimental  and 
theoretical  data  were  shown  to  correlate  well.  Time  histories  of  smooth 
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profile  slews  showed,  however,  that  the  shaping  did  not  work  as  well  for 
higher  modes  when  the  angle  of  the  slew  became  large  (>  1^).  Many 
“nonlinearities”  were  analyzed  to  be  the  cause  of  the  extra  vibration,  but  not 
one  was  found  to  be  a  large  factor.  It  was  concluded  that  the  large 
amplitudes  of  the  slew  cause  the  structure  to  no  longer  act  linearly.  The 
larger  slews  introduce  more  total  energy  into  the  system  and  cause  the 
structure  to  deflect  to  amplitudes  that  take  it  out  of  the  linear  approximation 
range.  We  then  ended  the  chapter  by  discussing  the  effect  of  higher  and 
lower  bandwidth  controllers  on  the  shaping  method.  It  was  determined  that 
the  shaped  inputs  require  more  precision  from  the  actuator,  and  therefore  the 
higher  bandwidth  controller  was  much  more  effective  in  implementing  the 
impulse  shaper.  The  lower  bandwidth  controller  was  not  only  ineffective,  but 
contributed  more  amplitude  to  the  residual  vibration  at  higher  frequencies. 

The  conclusions  drawn  in  Chapter  5  will  have  the  most  impact  on  impulse 
shaping  technology.  We  have  always  been  weary  of  how  impulse  shaping 
would  work  on  complicated  nonlinear  systems.  It  has  been  shown  here  how 
effective  it  can  be.  It  certainly  shows  an  improvement  over  unshaped  moves 
in  most  configurations  especially  in  the  lower  frequencies  But  as  can  be  seen 
from  higher  frequency  results,  it  doesn’t  meet  expectations  for  large  angle 
slews.  The  effectiveness  of  the  shaper  was  directly  correlated  to  the  slew 
angle.  We  saw  that  a  1°  slew  showed  similar  effectiveness  to  that  shown  by 
small  angle  frequency  response  experiments.  Larger  slews  exhibited  more 
shaped  residual  vibration.  Also,  a  correlation  was  found  between  the 
bandwidth  of  the  controller  and  the  effectiveness  of  shaping.  A  good  rule  of 
thumb  we  may  take  away  from  this  is  that  the  bandwidth  of  controller  should 
be  as  high  as  the  frequency  we’re  trying  to  cancel.  A  3Hz  bandwidth 
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controller  was  found  to  have  no  effect  on  trying  to  cancel  a  9Hz  mode,  while  a 
20Hz  controller  was  more  effective. 


Also,  we  have  seen  a  comprehensive  review  of  the  present  status  of 
impulse  shaping  theory.  The  tradeoffs  between  the  different  methods  are 
described,  and  suggestions  on  how  find  the  most  efficient  shaper  are  given. 
This  paper  will  hopefully  stand  as  a  resource  for  users  to  review  the  options 
available  in  implementing  impulse  shaping. 


6.2  Suggestions  for  MACE  Project 

I  regret  not  being  able  to  stay  on  the  MACE  project  to  see  it  to  completion. 
Seeing  years  of  work  come  together  would  have  been  a  gratifying  experience. 
I  am  comforted  by  the  fact  that  the  project  is  in  the  hands  of  some  very 
competent  and  talented  people  at  SERC.  I  would  like  to  leave  my 
impressions  and  suggestions  for  the  MACE  project,  however. 

As  discussed  in  chapter  3,  there  were  some  problems  with  the  gimbal 
hardware  and  the  AC-100  computer.  The  gimbal  had  structural  frequencies 
whose  amplitudes  were  high  enough  to  require  a  lowpass  filter  in  the  control 
loop.  Hopefully,  this  problem  will  be  solved  in  the  second  generation  gimbal 
design.  The  AC-100  computer,  although  a  very  user  friendly  environment, 
may  present  problems  with  processing  delay  when  dealing  with  multi-input, 
multi-output  systems.  For  example,  two  encoders  could  not  be  read  at  the 
same  time  at  the  current  lEHz  servo  rate.  Also,  some  of  the  more 
complicated  direct  solution  impiilse  sequences  could  not  be  implemented 
because  of  processing  delay.  A  new  processor  board  or  some  further 
investigation  by  Integrated  Systems  woiild  be  welcomed.  As  MACE  moves  on 
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to  higher  order  controllers,  the  stability  of  the  computer  will  become  more 
crucial.  Also,  the  original  specifications  for  the  slewing  size  of  the  gimbal 
have  not  been  met.  Presently,  the  gimbal  is  constrained  to  move  ±30**  by  its 
own  design.  The  suspension  system  constrains  that  movement  further  once 
the  gimbal  is  attached  to  the  bus  to  <  ±20**.  A  larger  slewing  angle  can  be 
more  useful  in  studying  the  non-linearity  of  the  system  and  the  gravity 
effects  of  a  cantilevered  payload.  Considerations  for  a  new  suspension  system 
or  a  lighter  payload  should  be  given. 

As  far  as  impulse  shaping,  much  can  be  learned  from  the  results  in 
chapter  5.  As  discussed,  a  high  bandwidth  controller  is  more  responsive  to 
the  needs  of  impulse  shaping,  and  hence  will  yield  more  of  the  vibration 
redudng  characteristics  of  the  shaper.  Since  th'  augh  dynamic  testing  of  the 
test  article  wiU  he  done  both  on  earth  and  in  sp^  ^e,  the  modal  firequendes  can 
be  identified  fairly  accurately.  Al^,  the  firequendes  do  not  seem  to  shift  very 
much  from  move  to  move,  at  least  in  the  range  of  motion  that  is  possible. 
This  may  suggest  a  shorter  sequence  in  exchange  for  lack  of  some  robustness. 
A  direct  solution  sequence  for  the  first  four  vertical  modes  plus  one  for  the 
two  horizontal  modes  around  4.7Hz  is  an  effective  and  relatively  short 
sequence  that  can  be  used.  Of  course,  these  firequendes  will  shift  as  different 
controllers  are  implemented.  If  the  structure  were  to  go  into  orbit  with  the 
current  controller,  which  it  probably  will  not,  this  is  the  impulse  sequence 
that  should  be  iised.  The  direct  solution  sequence  is  shown  below. 
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FIGURE  6.1:  A  five  mode  direct  soluti(xi  sequence  for  the  MA(I!E 
Develofunent  Model.  (1.75, 4.7, 52\  8.68, 9.14  Hz). 


0 

0.0313472 

0.0784381 

0.0665015 

0.155516 

0.0916246 

0.253105 

0.086278 

0.340187 

0.127523 

0.424599 

0.189451 

0.509011 

0.127523 

0.596093 

0.088278 

0.693681 

0.0916246 

0.770759 

0.0665015 

0.849198 

0.0313472 

TABLE  6.1:  Listing  of  impulse  times  and  amplitudes  fOT 
iiiq>ulse  sequence  in  figure  6.1. 
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Most  of  all,  experimental  data  needs  to  be  continually  gathered  on  various 
combinations  of  shapers  and  controllers.  Although  many  megabytes  of  data 
have  been  recorded,  there  are  so  many  variables  to  tune  that  it*s  hard  to 
imagine  that  a  full  optimization  has  been  achieved. 


6.3  Future  Work 

My  research  has  answered  many  of  the  questions  in  the  impulse  shaping 
field.  There  is  much  room  for  more  research  on  the  topic.  Some  suggestions 
are  presented  here 

As  shown  in  this  work,  the  controller  is  something  that  must  be  dealt  with 
in  implementing  shapers.  Further  investigation  into  this  topic  may  reveal  a 
method  of  integrating  the  design  of  the  controller  with  the  capabilities  of  the 
impulse  shaper.  For  example,  more  overshoot  firom  the  controller  may  be 
tolerated  because  of  the  vibration  reducing  characteristic  of  the  shaper.  A 
fiiamework  can  be  created  in  which  the  design  of  the  filter  and  the  controller 
can  be  simultaneously  addressed.  An  even  further  extension  of  this  topic 
may  be  to  incorporate  the  shaper  inside  the  control  loop.  Currently  this  is  not 
possible  due  to  the  phase  lag  inherent  in  an  impulse  shaper.  A  closed  loop 
implementation  will  work  not  only  to  cancel  vibration  caused  by  the  move, 
but  also  vibration  caused  by  disturbances  in  the  system. 

Impulse  shapers  work  to  cancel  residual  vibration  but  do  not  affect  the 
vibration  during  the  move.  In  tracking  applications  in  space,  the  vibration 
during  the  move  is  equally  important.  A  new  formulation  of  the  shapers  that 
addresses  this  issue  could  be  helpful  to  the  MACE  project. 
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An  investigation  into  a  non-linear  extension  of  impulse  shaping  would  be 
a  rewarding  research  area.  The  results  from  chapter  5  showed  that  non- 
linearities  limited  the  effectiveness  of  the  shaper  for  higher  firequency  modes. 
Although  some  suggestions  were  made  to  help  the  problem,  the  non- 
linearities  were  not  directly  addressed.  A  non-linear  formulation  of  the 
vibration  constraint  equations  or  an  extension  of  the  existing  theory  to 
approximate  non-linear  behavior  would  be  a  worthy  study. 

Further  investigation  of  the  fonnulation  of  direct  solution  multiple  mode 
sequences  is  needed.  The  complexity  and  sometimes  inability  of  the  method 
to  generate  solutions  makes  it  difficult  to  use  for  some  applications.  A  closed 
form  solution  of  this  method  would  prove  very  beneficial. 

Finally,  more  comparison  studies  need  to  be  conducted  to  illustrate  the 
advantages  of  impulse  shaping  over  other  types  of  input  shaping.  This 
examination  could  lead  to  developments  of  shaping  algorithms  that  utilize 
many  types  of  input  shapers.  A  shaper  that  applies  impulse  shapers  for  low 
modes  and  digital  notch  filters  for  high  modes  is  an  example.  It  is  clear  that 
the  good  attributes  of  the  different  types  of  input  shaping  can  be  brought 
together  to  maximize  the  efficiency. 
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